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We consider Hamiltonian systems which can be described both classically and quantum mechan- 
ically. Trace formulas establish links between the energy spectra of the quantum description, and 
the spectrum of actions of periodic orbits in the Newtonian description. This is the duality which 
we investigate in the present paper. The duality holds for chaotic as well as for integrable systems. 

Billiard systems are a very convenient paradigms and we use them for most of our discussions. 
However, we also show how to transcribe the results to general Hamiltonian systems. In billiards, 
it is natural to think of the quantal spectrum (eigenvalues of the Helmholtz equation) and the 
classical spectrum (lengths of periodic orbits) as two manifestations of the properties of the billiard 
boundary. The trace formula express this link since it can be thought of as a Fourier transform 
relation between the classical and the quantum spectral densities. It follows that the two point 
statistics of the quantal spectrum is related to the two point statistics of the classical spectrum 
via a double Fourier transform. The universal correlations of the quantal spectrum are well known, 
consequently one can deduce the classical universal correlations. In particular, an explicit expression 
for the scale of the classical correlations is derived and interpreted. This allows a further extension 
of the formalism to the case of complex billiard systems, and in particular to the most interesting 
case of diffusive system. The effects of symmetry and symmetry-breaking are also discussed. 

The concept of classical correlations allows a better understanding of the so-called diagonal ap- 
proximation and its breakdown. It also paves the way towards a semiclassical theory that is capable 
of global description of spectral statistics beyond the breaktime. An illustrative application is the 
derivation of the disorder-limited breaktime in case of a disordered chain, thus obtaining a semi- 
classical theory for localization. We also discuss other applications such as the two-cell systems, 
periodic chains and localization theory in more than one dimension. 

A numerical study of classical correlations in the case of the 3D Sinai billiard is presented. Here 
it is possible to test some assumptions and conjectures that underly our formulation. In particular 
we gain a direct understanding of specific statistical properties of the classical spectrum, as well as 
their semiclassical manifestation in the quantal spectrum. 

We also analyze the spectral duality for integrable systems, and show that the Poissonian statistics 
of both the classical and the quantum spectra can be traced to the same origin. 



I. INTRODUCTION 



Trace formulas manifest an intimate link between two seemingly unrelated spectra which pertain to the quantum 
and the classical descriptions of a given Hamiltonian system. The quantum spectrum of eigenenergies is described in 
terms of the spectral density 

oo 

dg^{E;n) = ^^[i;-£;„(n)] (i.i) 

n=l 

and it depends parametrically on Planck's constant. The classical spectrum of actions characterizes the set of periodic 
manifolds at an energy E. These are isolated periodic orbits (POs) in case of a chaotic system or periodic (rational) 
tori in case of integrable system. The actions are expressed as contour integrals in the 2d dimensional phase space 

„ d 

Sa^iE) = j> ^p.dg, (1.2) 
•^C"" i=l 

where Ca is a closed trajectory on the periodic manifold (PO in the case of a chaotic system). Note that the periodic 
manifolds, as well as their actions, depend parametrically on the chosen energy E. One can define a density of classical 
actions 

d,i{S-E) = ;^A«^[5-5„(£;)] (1.3) 
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and with an appropriate choice of the coefficients Ac the trace formulas are expressed as 



,iE;h) ~ (27r?i) 



{27rh)-f j e^/''d,i^osc{S-E) dS 



(1.4) 



where dqm,osc = dqm — dqm and dqm{E; h) is is the smoothed spectral density. The power / depends on the dimen- 
sionality of the classical periodic manifold: it is 1 for isolated POs in chaotic systems, and (d + l)/2 for periodic tori 
in integrable systems. The symbol ^ indicates that the right hand side provides the leading term in an asymptotic 
expansion in fi of the left hand side. This is the semiclassical trace formula (SCTF) which is the main object of 
our discussion. Explicit expressions for the semiclassical Aa coefficients for chaotic and for integrable systems were 
derived by Gutzwiller Q and Berry and Tabor respectively. It should be noted, that there are few cases where 
the symbol can be replaced by an equality sign. One example is the Selberg trace formula for the modular domains 
in the hyperbolic plane. A second example is the rectangular billiards in any dimension. The higher terms in the Ti 
expansion are known in some cases [ pT] and are introduced by adding to A^ a power s eries in fi. For a general billiard 
there is a theorem by Andersson and Melrose H that states that exist A^ for which (|l.4|) is an exact equality. 



It should be emphasized that both sides of the SCTF (1.4) are distributions rather than functions. Therefore one 
cannot use it directly but only after it is integrated over an appropriate test or window function. Once this is done 
one can use the SCTF to study the quantal spectral statistics of a finite spectral interval in terms of the properties 
of the classical dynamics which are embedded in dd ||l^,|l^,^ . Consider the two point spectral form factor, which is 
calculated for the spectral window w{E — E) whose effective width is Ai?. For t > 



K{t-E) = 



dqmAE 



iE„t/h 



(1.5) 



The inner summation over n can be replaced by an integration over E with the density dqm{E). Substituting the 
SCTF for dqm{E) and assuming that Ai? is sufficiently small on the classical scale, but large enough on the quantal 
scale, one obtains: 



K{t-E) 



-f 



dqmAE 



^.^(t-r„)A„e'^°(^)/'' 



(1.6) 



Here w is the Fourier transform of the energy window function, and it pick s out of the PO sum only those POs whose 
period approach t to within h/ AE. Note that the sum over POs in (1^) absolutely converges for a suitably chosen 
window. The absolute square of the sum in (1.6) is composed of two kinds of contributions. The contribution of the 
diagonal terms is 



KD{t,E) 



{2TThy 



dqmAE ^ 



^'w{t-T^)f\A^f 



(1.7) 



and the contribution of the non diagonal terms is 

K^o{t- E) = ^^-""^l I V w{t - T^)w{t - Tp) A^Al e^is^-s.]/"' 



(1.8) 



Berry ||T^ was the first to observe that for times much shorter than the Heisenberg time tniE) = 2-Khdqm the 
diagonal contribution dominates. Using the Hannay and Ozorio de Almeida sum-rule JT^ he was able to reproduce 
the expression of K{t) for t < tn as derived from random matrix theory (RMT) for generic chaotic systems, and the 
expression derived for a Poissonian spectrum for generic integrable systems. This approach has been generalized later 
by observing that the diagonal sum is r elate d to the classical probability to return fl^ . 
The contribution of the non diagonal sum ( |1.§| ) is by no means small for t > tn- As a matter of fact it should contain 
a term which (for chaotic systems) cancels the monotonically increasing Kjj cx t so that the correct asymptotic behavior 
K(t) 1 is achieved for t > tn- This asymptotic behavior refiects the fact that the density dqm{E) corresponds 
to a discrete spectrum. Note that if dqm{E) were smooth, then K{t) would asymptotically vanish. The asymptotic 
non-zero value is therefore a hallmark of the discreteness. Argaman et al. have related K^jjit) to the two 
point correlations of the classical spectrum. The distribution of the classical actions differences [Sa — Sp] , with the 
restriction tha t th e periods are confined to the vicinity of the time t, has been determined by taking an inverse Fourier 
transform of ( |1.8D in the variable 1/h. For chaotic systems in 2D, the classical two-point correlations scales into a 



universal function, provided the action differences are expressed in units of j 



on I 



t dE\E' 



where fl{E) is the volume of the 
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phase space with energy < E. This prediction was checked numerically for a few chaotic systems, and its relation to 
the Hardy-Littlewood conjecture on pair correlations of primes was also discussed . 

The important new element in was that a hitherto unnoticed correlations between classical actions were derived 
from the SCTF, and from the phenomenological observation that the spectral statistics of quantum chaotic systems 
reproduce the results of RMT. (See also the discussion by Dahlqvist ^^). It should be emphasized that it is crucial to 
assume that the SCTF indeed reproduces discrete energy levels ('convergence' in the sense of distributions), else the 
argument for the existence of 'classical correlations' loses its validity. Thus, there is an inherent connection between 
the quest for classical correlations, and the wonder concerning to the mathematical soundness ('convergence') of the 
SCTF. This point has been raised by Gutzwiller Q who introduced the term 'third entropy' in order to refer to the 
statistical properties of the actions. 

The conjecture of classical action-correlations leads naturally to the introduction of some fundamental questions: 

• What is the dynamical origin of the correlations between classical actions? In other words, how can one derive 
these correlations from purely classical considerations, without recourse to the SCTF and quantum mechanics ? 

• The classical density consists of S functions which are supported by the actions of the periodic manifolds, and 
weighted by complex coefficients whose magnitude as well as phase differ from orbit to orbit. How much of the 
expected correlations reflect the distribution of actions, and what is the role (if any) of the correlations amongst 
the coefficients themselves or amongst the actions and the coefficients? 

• Periodic orbits proliferate exponentially with their periods. Excluding non generic degeneracies, this implies an 
exponentially large density of actions of POs. Previous numerical checks |13| have convincingly shown that on 
the scale of the (exponentially small) nearest actions separation, the action spectrum is Poissonian. How is this 
compatible with action correlations? 

• What is the classical phase-space significance of the classical correlation scale (geometrical interpretation in case 
of billiards)? What is the way to generalize the scaling relation for action correlations beyond 2D ? What are 
the modifications that are required in case of complex (e.g. diffusive) systems ? 

Even though the non diagonal contribution to the form factor attracted much attention in the past few years, all the 
studies so far avoided the direct confrontation with action correlations and the questions listed above. Khmelnitskii and 
coworkers ||2^ as well as Altshuler et al computed the form factor for disordered (chaotic) systems by employing 
techniques other than the semiclassical theory. Bogomolny and Keating ||2^] effectively replaced the contributions 
of POs with t > by composite POs constructed in a special (synthetic) way. Miller has deduced action 

correlations of composite POs from unitarity, and derived from them the classical correlations of primitive POs. 

The purpose of the present paper is to pick up again the subject of action correlations from the point where it 
was left by jiot . We shall introduce an explicit formulation of the quantal- classical duality concept, and will extend 
it to the general case of either integrable or chaotic systems with d degrees of freedom. We shall refer to generic 
as well as to non-generic circumstances. At a later stage we shall discuss the modifications that are required in 
the application to complex (e.g. diffusive) systems. The formulation of quantal-classical duality and the associated 
heuristic understanding of classical correlations thus paves the way towards a global semiclassical theory of spectral 
statistics that goes beyond the diagonal approximation. We use from now on the term 'classical correlations' or 
'classical PO-correlations' rather than the term 'action correlations'. Sometimes, the latter term is misleading. Some 
of the results were already reported by D. Cohen ||2l[| . They are discussed and further developed in the present paper. 

Some remarks on the method of presentation are in order. Most of the discussion will be carried out for billiard 
systems. Billiards are particularly convenient since they are scaling systems. Because of this property, the SCTF can 
be simplified by considering the quantum wavenumbers instead of the energies, and in the classical description the 
length of the PO replaces the action and the time, (since one can always choose m = v = 1). Once the theory for 
billiards is presented, we shall show how to transcribe it for any Hamiltonian system. 

In section II the general formulation of the quantal-classical duality is presented. The SCTF is used in order to 
derive the relation between the two point statistics of the quantal and the classical spectra. In particular, the universal 
classical correlation scale A is identified. This section contains a detailed discussion of the statistical procedure and 
its relation to the energy-time plane. We discuss the crossover from the 'classical regime' to the 'quantal regime' 
and the proper definition of the breaktime follows naturally. The applicability of the semiclassical theory beyond the 
(irrelevant) Ehrcnfest 'log'-time is clarified. The role of non-universal features is discussed as well. 

In section III we give an interpretation of the classical correlation scale, and introduce a statistical model for 
rigidity in terms of 'families'. Then we use various arguments and the notions of 'classification' and 'similarity' in 
order to reveal the actual structure of the classical spectrum. Finally we make some observations concerning the 
origin of classical correlations. 
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Some of our observations concerning classical correlations are further supported by the numerical study of section 
VI. There we have carried out a detailed analysis of the quantum and classical spectra of the Sinai billiard in 3D, for 
which we possess substantial and immaculate data bases |^^. This analysis provides for the first time a clue to the 
basic problem - the derivation of the classical PO-correlations based on a purely classical argument rather than on 
the semiclassical duality. 

Section IV extends the formulation of section II to the general case of either chaotic or integrable Hamiltonian 
systems. For generic integrable system we discuss the dual Poissonian nature of both the quantal and the classical 
spectra. (See also |^-|l0|). The general strategy is further illustrated by referring to specific examples, namely, the 
Riemann zeros and the 3D torus. The cubic 3D torus is particularly interesting due to its non-genericity and is 
analyzed in detailed. 

Besides the mathematically oriented questions that concern the existence of classical correlations, there is also 
the practical hope to extend the semiclassical theory beyond the diagonal approximation. Classical PO-correlations 
become most significant in the analysis of complex systems since they lead to relatively short breaktime scales that 
are not related to the universal Heisenberg time. In section V we make use of the physical insight gained in the 
discussion of classical PO-correlations in order to formulate the modifications that are required in order to go beyond 
the diagonal approximation. A specific example is the quasi-lD classically-diffusive chain. The quantal form factor 
should correspond either to Anderson-localized eigenstates or to band structure, depending on whether the chain is 
disordered or periodic. These genuine quantal effects are due to interference, hence the non-diagonal terms play a 
crucial role. One may question the feasibility of having PO-theory for the spectral properties in such circumstances. 
Thus, our interest in complex systems is twofold. On the one hand we develop a semiclassical theory which is not 
limited by the diagonal approximation. On the other hand we test the limits of the semiclassical approach. We 
use our present level of understanding of the SCTF and of classical PO-correlations in order to demonstrate that the 
SCTF should be effective also in the study of complex systems. The emergence of a non-universal volume-independent 
breaktime in case of a disordered systems is an obvious consequence of the formulation. 



II. GENERAL THEORY 



The semiclassical trace formula (SCTF) relates the quantal spectrum of energy levels to the classical spectrum of 
periodic orbits (POs). As explained in the introduction the SCTF can be utilized in order to derive a semiclassical 
expression for the quantal form factor. The purpose of this section, following [ pl| , is to introduce an explicit formulation 
of the quantal-classical duality concept. The reader should notice that a clear and explicit formulation of the quantal- 
classical duality requires modifications of some common notations. In particular, we find it convenient neither to 
unfold the quantal spectrum nor to scale the quantal form factor. It is also most convenient to adhere to the standard 
Fourier transform conventions. The new formulation, which will be discussed thoroughly in the following section, 
promotes the understanding of the classical correlation scale, illuminates the crossover from the 'classical' to the 
'quantal' regime, and sheds new light on the significance of the statistical procedure. 



A. The Semiclassical Trace Formula 



The quantal spectrum {fc„} consists of real positive eigen-wavenumbers that are defined by the Helmholtz equation 
(V^ -I- k'^)ip = with the appropriate boundary conditions. The average (smoothed) density is found via Weyl law 
whose leading order term is 



k'^ 



Sk 



Cd = 



Z TT - 



1 

-^r(-) 



(2.1) 



In the above fl is the volume of the billiard, and d is the dimensionality. The statistical averaging (...) implies smoothing 
using a window of width 6k. Further discussion of the statistical procedure will appear at later subsections. It is 
natural to define a density Pqm{k) that corresponds to the quantal spectrum, namely 



Pqm{k) = 2TT[dqm{k) 



rfqm(fc)] = 27r(5(fc-fcn) 



(2.2) 



The subscript osc implies that the smooth component of the density is being subtracted. This is done in order to 
simplify the subsequent formulation. Furthermore, in order to facilitate the application of Fourier transform (FT) 
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conventions we define p{—k) = p{k) for /c < 0. The factor 27r has been incorporated in the above definition for the 
same reason. 

The classical spectrum {Lj} of periodic orbits (POs) is defined as the set of all the primitive lengths Lp and their 
repetitions 0<r such that Lj = rLp. For a simple chaotic billiard, due to ergodicity, the so called 'Hannay and Ozorio 
de Almeida sum rule' implies 



Ku{L) ^ IY,\A,\H{L-L,) 



L 



(2.3) 



SL 



where a smoothing window of width 5L is being used. The generalization of this sum rule to complex systems is 
discussed later in section V-B. The weighting factors are \Aj\'^ = Lp/\ det{Mp—I)\ where M is the monodromy matrix. 
These weighting factors decay exponentially with L, namely \Aj\'^ ~ exp(— crL), where a is the average Lyapunov 



exponent. Hence (using ( p.3[ )) one deduces that the actual smoothed density dci{L) of POs grows exponentially as 
exp{aL)/L. It is natural to define a weighted density Pci{L) that corresponds to the classical spectrum, namely 



A.- = p 

y^|det(M;-/)| 



(2.4) 



rVp which 



The instability amplitudes Aj are endowed with a phase factor that contains the effective Maslov index Vj 
includes repetitions and the boundary contribution due to reflections. As in the analogous definition of the quantal 
density, we adhere to standard FT conventions: The classical density satisfies the symmetry relation p{—L) — p*{L), 
and the sign in front of the Maslov index is positive. A negative sign would be incorporated if time domain FT 
conventions (—iutt) rather than space domain FT conventions {+ikx) were used. The former are used in the general 
derivation of Gutzwiller . 

The semiclassical trace formula (SCTF) relates the quantal density Pq„i{k) to the classical density pci(L). The 
relation is simply 



Pqmik) = TT Pcl{L) 



(2.5) 



Disregarding the smooth part, the SCTF (2^) states that ^„ 27r(5(fc— fc„) = ^^^^ exp(— ifcLj). Its Fourier trans- 
formed version is ^„exp(ifc„L) = J^j AjS{L~Lj). However, if the two sides of the trace formula were multiplied by 
a window function then the Fourier transformed version would be (See remark [Q), 



E 



w{kn—k) exp(iA;„i) 



^kL ^ w(^Lj-L)Aj expi^ikLj) 



(2.6) 



where w is the Fourier transformed window. This version of the trace formula is manifestly symmetric and is also 
more convenient for numerical studies. From a mathematical point of view the latter version of the SCTF is superior 
due to its convergence property. 



B. Duality of Two Points Statistics 



In order to study the statistical properties of either the quantal or the classical spectrum it is natural to define 
the corresponding two point correlation functions. It follows from the SCTF that the two-point statistics of the 
quantal spectrum is related to the two-point statistics of the classical spectrum via a double Fourier transform, 
namely {pqm{k)pqm{k')) = TTFT {pci{L)pci{L')). Without the averaging operation (...), the latter relation is 
mathematically trivial and useless. In the sequel, we shall argue that the equality holds non-trivially as a statistical 
relation, meaning that the averages on both sides of the equality can be carried out independently. We postpone for 
a moment both the discussion of this point, and the precise definition of the averaging procedure. We turn first to 
introduce some formal definitions. 

The two point correlation function of the quantal spectrum is defined as follows 

(2.7) 

^ - (2W,„(fc))2 (2.8) 

(2.9) 



Rqm{k,e) = \Pqm(k-^j Pqm (k+^jj 
= ( E H^~ikn^-kn)) S(k 

= 27rKQ{k)-[d{e)-b{e;k)] . 
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It is a function of e and it depends parametrically on k. Sometimes it is convenient to notationally suppress the 
parametric dependence. We shall frequently use the brief notation Rqm(e)- The smoothed density squared has 
been subtracted, and therefore Rqm{e) approaches zero asymptotically for e ^ oo. The i5(e) term is due to the self 
correlations (n = m) of the levels. The function 6(e) can be interpreted as the 'probability' density for finding a 
vacancy at a range e from some reference level. It should have normalization equal to one since the quantal spectrum 
is characterized by a finite rigidity scale, still it may possess negative parts if there is a clustering effect. The spectral 
form factor Kqm{k, L) is the FT of Rqrn{k, e) in the variable e^L. 

The two point correlation function of the classical spectrum is defined analogously as follows 



Rci{x,L) 



Pel 



Kd{L)-[5{x)~p{x-L)] 



(2.10) 
(2.11) 
(2.12) 



It is a function of x and it depends parametrically on L. Sometimes it is convenient to notationally suppress the 
parametric dependence on L. The smoothed density squared has subtracted, and therefore Rci{x) approaches zero 
asymptotically for x — *■ oo. The S{x) term is due to the self correlations of POs. The function p{x) can be interpreted 
as the 'probability' density for finding a vacancy at a range x from some reference PO. It should have normalization 
equal to one if the the spectrum is characterized by a finite rigidity scale. We shall argue later that this is indeed the 
case. Still, it may possess negative parts if there is a clustering effect. 

With the above definitions the relation between the the quantal two points statistics and the classical two points 
statistics can be written in the following form 



Rqm{k,e) = KTTT Rci{x,L) 
where the first !FT corresponds to x 



(2.13) 



k and the second JH' corresponds to L ■ 



Alternatively 



qUk,L) = Kci{k,L) . (2.14) 

Thus, on semiclassical grounds, if the function K(k, L) is viewed as a function of L then it equals the quantal form 
factor, while if it is viewed as a function of k it equals the classical form factor. It also follows that if the quantal 
spectrum is characterized by non-trivial correlations, then also the classical spectrum should be characterized by some 
non-trivial correlations. 

We turn now to define the meaning of the statistical averaging. It is most convenient to do it in the (fc, L) plane. 
Namely, 



{f{k,L))sk5L = 



1 



SkSL 



f{k',L') dk'dL' 



(2.15) 



In the {k,e) plane the SL averaging implies multiplying Rqmif) with a cutoff function whose width is 271 /5L. Analo- 
gously, in the {x,L) plane the Sk averaging implies multiplying Rci{x) with a cutoff function whose width is 2T:/dk. 
The semiclassical relation holds trivially provided the same window parameters are used in both sides of the equation, 
namely 



Kqm{k, L,Sk,SL) = Kciik, L,5k,SL) 
If the further condition 2n ^ SkSL is satisfied, then this relation can be cast into the form 



(2.16) 



y^exp(ifcnL) 



Aj exp{ikLj 



-<Lj<L + 2 



(2.17) 



The condition 27r <^ SkSL guarantees that the two points statistics is not sensitive to the exact value of the window 
parameters. It follo ws that we may choose the window parameters that appear in both sides of ( 2.16| ) independently. 
In this sense ( 2.14 ) becomes a statistically meaningful relation. Further discussion of the {k, L,Sk,6L) space will 
appear at later subsections. 

The summations in the above statistical relation are the same which occur in the calculation of the SCTF. The 
number of terms in these summation s is finite as in (2.6), due to the windows which are used. Equation ( ^.17 ) is in 
some sense the 'squared' version of (2.6), the main difference being that in the latter case 6L = on the left hand 
side, while dk = on the right hand side. Furthermore, in (2.6) the L-window shou ld be related to the fc-window by 
a Fourier transform. Later we shall discuss one more aspect of the relation of ( ^.17 ) to the SCTF. 
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C. Universal correlations 



The quantal spectrum is characterized by a universal correlation scale A(fc). For a simple ballistic billiard it simply 
coincides with the average level spacing. Random Matrix Theory (RMT) further determines the appropriate universal 
scaling function, namely 



6(e;fc) = 



1 



(2.18) 



A(fc) \A{k)^ 

From this scaling relation one derives the corresponding scaling relation of the quantal form factor, namely 

Kq^ik,L) = LH{k)- ko{L/LHm , (2.19) 

where i//(fc) = 27r/A(fc) = KQ{k) is the Heisenberg time. For the Gaussian Unitary Ensemble (GUE) there is a 
particularly simple result Kq{t) — min(T, 1). If there is a time reversal symmetry one should use the result of the 
Gaussian Orthogonal Ensemble (GOE). 

Using analogous strategy, let us assume that the classical spectrum is characterized by some universal correlation 
scale A(L) such that 

pix;L) = ^p(^) , (2.20) 



where p{s) is some scaling function. The justification for assuming the above scaling behavior will be discussed in 
subsection F. For simplicity of presentation it has been assumed that the classical spectrum has no degeneracies. We 
note however that if time reversal symmetry is taken into account, where for the generic non self tracing orbits the 
degeneracy is g = 2, then the scaling relation should be modified by writing p(x) — ~{g~l)S{x) + [g / X)p{x / \) so that 
Rci{x) = g[S{x) — {1/X)p{x/X)]. One may say that this degeneracy constitutes a trivial type of classical correlations. 

Taking the Fourier transform of Rqm{i) we have deduced that Kq„i{L) should have a crossover at the Heisenberg 
time L = 2tt/ A{k). Similarly, taking the Fourier transform oi Rci{x) it is obvious that Kci{k) should have a crossover at 
k = 2-k/\{L). Alternatively, one may say that the crossover should occur when \{L) equals the De-Broglie wavelength 
2'K/k. However, due to the semiclassical relation ( 2.14 ) the latter condition should coincide with the former condition 
for having a crossover at the Heisenberg time. It follows that K{k, L) has two regimes in the (fc, L) plane which are 
separated by a crossover line whose equation is 

27r , 27r , ^ 

L=ir7TT ^ k=—— . 2.21 
A(fc) A(L) ^ ' 

Assuming that A(fc) and A(L) are monotonic decreasing functions, it follows that the {k,L) plane is divided into a 
classical low-L-large-fc regime (C), and a quantal large-L-low-fc regim e (Q). These are illustrated in the left drawing 



of Fig.0. Using Weyl law 27r/A(fc) — Cd^k*^ ^, it follows from ( 2.21 ) that for a simple ballistic billiard the classical 
correlation scale should be 

ML,^) = 27rfQ-j (2.22) 

Note that this prediction is independent of any detailed RMT result. For a 2 D syst em, it coincides with the scaling 
relation of jl^ that we have mentioned in the introduction. The validity of ( 2.22|) follows from the validity of the 
SGTF. The subscript indicates that the universal (rather than some non-universal) correlation scale were considered 
here. 
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FIG. 1. Regions in the (k, L, Sk, 5L) space. The classical (C) and the quantal (Q) regimes, as well as the classical-diagonal (CD), 
classical-statistical (CS), quantal-diagonal (QD) and quantal-statistical (QS) sub-regimes are indicated. The mean chord is C. Individual 
POs and energy levels are represented in the left drawing by small bars. The dashed line is the SkSL = 2it border for statistical stability. 
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FIG. 2. Schematic plots of -Rqm(e) and Kqm{L) for k = const, and of K^i{k) and Rci{k) for L = const. The dashed Unes reflect the 
quantal diagonal approximation, and the dotted lines reflects the classical diagonal approximation. Kqm{L) is obtained from Rqm{^) via 
FT, Ki-i = Kqm from semiclassical considerations, and R^-i is obtained from K^-i via an inverse FT. Both the quantal (Q) domain and the 
classical (C) domain are indicated in the plots of the form factor. The origin of the vertical axis in the plots of the Rqm and Rd is shifted 
upwards since we subtract the smooth background component. 



D. The Diagonal approximation 

We turn now to discuss the so called "diagonal approximation" for the form factor. By diagonal approximation it 
is meant that the non trivial correlations of either the quantal or the classical spectrum are ignored. In the absence 
of degeneracies only self-correlations are taken into account. For the quantal form factor, the diagonal approximation 
is 

K qm{k , L) K, K Q{k) [27r/A(fc)<Ci] Diagonal approximation, deep quantal regime . (2.23) 
The analogous diagonal approximation that refers to the classical form factor is 

Kci{k, L) fi: K]j{L) [27r/A(i)^fc] Diagonal approximation, deep classical regime . (2.24) 
If the classical spectrum has an effective degeneracy g due to some symmetry then one should write Kci{k,L) 



gKu{L). Thus, upon using the semiclassical relation (2.14), the asymptotic behavior of the function K{k,L) is 



determined both in the deep classical regime (C) and in the deep quantal regime (Q). This is illustrated schematically 
in Fig.| 

The above approximations may become an exact identities if appropriate windows are used. In the deep quantal 
regime one may use a window such that 27r/A ^ 6L. It is equivalent to multiplication of Rqm{i) by a cutoff window 
which is much narrower than the correlation scale A. It implies that the non-trivial correlations that are represented 
by 6(e) are washed out, and therefore Kqm{k, L) becomes identical with the one-point function KQ{k). The sub-regime 
in (fc, L, 6k, SL) space where the latter observation applies will be called the Quantal-Diagonal (QD) regime. Similarly, 
we may define a Classical-Diagonal (CD) sub-regime. The latter occurs at the deep classical regime if 27r/A <C Sk. In 
this latter sub- regime one may argue that Kci{k, L) becomes identical with Ko{L). The different sub-regimes in the 
{k, L, Sk, SL) space are demonstrated in Fig.|l|. 



In both the QD and the CD regimes the semiclassical relation (2.14) obviously holds. However, it looses some of its 
statistical significance since it does no longer represent a two-way statistical relation. Rather it reduces to a 'one-way' 
statistical relation where one side of the equation is a one-point entity. An extreme case would be such that either 
Sk is smaller than the average level-spacing Afc, or SL is smaller than the average length-spacing AL. The latter 
conditions are sometimes assumed while analyzing non-universal long range correlations in the quantal spectrum. 

It is interesting to note the behavior of the classical form factor Kci{k,L) in the quantal diagonal (QD) regime. 
The summation in the right hand side of ( 2.17] ) is the same which occurs in the calculation of the SCTF. It will yield 



delta functions of width 211 /SL. These delta functions are resolved (at least in the statistical sense). After squaring, 
it should give a semiclassical estimate for the one-point density of the quantal spectrum. Analogous behavior is found 
for the quantal form factor in the CD regime. 



E. Breaktime Concept 

The semiclassical description of quantum-dynamics consists of several time regimes. In this subsection we wish to 
make a conceptual distinction between time scales which are associated with the breakdown of the stationary phase 
approximation, and those which are associated with interference effects. The latter are important while discussing 
the diagonal approximation for the form factor. 

Important Note: From now on we shall frequently translate lengths into times by using L = vt where v is the 
velocity. Similarly one may translate Ak into conventional energy units by using AE = vAk. Formulas and relations 
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such as ^i: I ^ < L are easily transformed into the corresponding time domain version 27r/A < t provided the proper 
'units' are used for A, or any other relevant object. As far as the semiclassical limit is concerned k plays the role of 1/h. 
Using different terminology one may say that classical mechanics is concerned with lengths of rays, while quantum 
mechanics associates with these rays a wavenumber k. The mass and the velocity of the particle by themselves are 
insignificant physically. The velocity v is used merely for converting units of length into units of time and therefore 
can be set equal to 1. 

The correspondence principle implies that quantal evolution should follow the classical one on short time scales. 
Deviation due to breakdown of the stationary phase approximation are expected after the time tsd- It has been argued 
that tsci ~ k^^^. Further deviations from the leading order semiclassical expansion due to e.g. diffraction effects 
are discussed in . One should be careful not to confuse these deviations which are associated with the accuracy of 
the stationary phase approximation with the following discussion of the breaktime concept. We assume in this paper 
that the leading order semiclassical formalism constitutes a qualitatively good approximation also for tsd < t in spite 
of these deviations. 

Interference effects lead to further deviations from the classical behavior. Well isolated classical paths that are 
involved in common semiclassical calculations may give rise to either constructive or destructive interference effects. 
In particular, interference contribution may be significant in the semiclassical computation of the form factor (right side 
of ( ^.17] )). The diagonal elements of the sum J^ij ^i^j 'Sxp{ik{Lj—Li)) represent the classical contribution due to self- 
correlations of POs, while the off diagonal part of the sum constitute the interference effect. The actual contribution 
of interference is determined by the statistical properties of the classical spectrum. If the classical spectrum is of 
Poisson type, then the interference contribution is self-averaged to zero. We shall use the term breaktime (t*) in order 
to denote the relevant time scale for the manifestation of such interference effect. 

If only the universal classical correlations are considered than t* should be identified with the Heisenberg time scale 
tn- Recall that the Heisenberg time is 27r/A(/c), hence tjy ~ k'^~^ which is semiclassically much larger than tsd- For 
t tn the diagonal approximation ( 2.24| ) should hold. Deviations from the diagonal approximation may be apparent 
in the classical (C) regime t < tn, and depend on the actual functional form of p{x). For GUE, for some unknown 
reason, the diagonal approximation actually holds over the whole classical (C) regime. In the quantal (Q) regime 
tn < t, the diagonal approximation fails completely. On time scales tn <^ t the recurrent quasiperiodic nature of the 
dynamics is revealed. 

If one confused the statistical scale Aq with the classical spacing AL Lexp(— crL), then one would obtain a 
false condition 2n <^ kAL for the validity of the diagonal approximation. Such confusion would lead to the wrong 
conclusion that there should be a breaktime at the 'log' time ts ^ In(fc), also known as the Ehrenfest time. The 
time scale ig, over which classical POs proliferates on the uncertainty scale, is semiclassically much shorter than 
both tsd and tH- The false condition 2tt <^ kAL is neither necessary nor sufficient for the validity of the diagonal 
approximation. This point is further discussed in the introduction for section III, and subsection III-F in particular. 
Summarizing, the message is that the 'log' time ts has no physical significance as far as the spectral form-factor is 
concerned. It is neither related to the breakdown of the stationary phase approximation, nor to the interference that 
leads to the breakdown of the diagonal approximation. 

In case of complex (non-generic) billiard system, the classical spectrum may be characterized by a non-universal 
shorter correlation scale A_. Such occurrences will be discussed in later sections. As a result there may appear a 
distinct breaktime t* which is shorter than tn- The crossover to quasiperiodic behavior may occur either at t* or still 
at tn- In the former case tn may loss its physical significance, as in the theory of quantum localization pi[|. 



F. Beyond the diagonal approximation 



In order to describe the departure from the diagonal approximation we shall define a correlation factor C(fc, L) via 
the relation 



Kdik,L) = C{k,L)-KD{L) 



Hence 



C{k,L) 



p{x;L) e 



-ikx 



dx 



(2.25) 
(2.26) 



If the classical spectrum has no degeneracies then the correlation factor should equal one in the deep classical regime 
L <^ 27r/A(fc). There, the classical diagonal approximation is expecte d to h old. The correlation factor should go to zero 
in the deep quantal regime 27r/A(fc) <^ L. This si niply follows from ( 2.25D using the fact that Kqm{k, L) KQ{k) is 
finite, while K£i{L) ^ L oo. It follows also (via (2.26)) that p{x) should have the normalization one, implying finite 
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rigidit y sca le. Thus it is natural to introduce a scaling function p{s) with scaling behavior that reflects normalization 
as in ( 2.20| ), and with scaling constant A = Xo{L,i}). It follows that there is a related scaling function such that 



C{k,L) = Co{kXo{L,fl)). More generally, if the the classical spectrum has degeneracy g than the proper scaling 
relation is C(fc, L) = gCo{kXo{L, fl)). This will be further discussed in section V-A. 

The considerations above actually determines the functional form of the scaling function Co(k) both in the deep 
classical regime I ^k a nd also in the deep quantal regime k^I. The interpolation requires further information. Using 



the RMT result (2.1£) it follows that the classical scaling function Co(k) is related to the quantal RMT scaling 




function Ko{t). For GUE statistics the scaling relation is C{k,L) = Co(fcAo(i, ri)), and the scaling function is 

for K^27r [corresponds to discrete quantal spectrum] 

for K<27r [RMT based interpolation - GUE ensemble] „_n 
for 2tt<k [RMT based interpolation - GUE ensemble] ' 
for 27r^K [based on the diagonal approximation] 

where the subscript has been added in order to emphasize that universal correlations are discussed. For GOE 
statistics the scaling relation is C{k,L) ~ gCa{kXg{L,n)), with Xg{L) = Xo{L,il/g) and the scaling function is 

(^)'^^^ f'-"' '^^Stt [corresponds to discrete quantal spectrum] 

^^^^^ _ . 2(^)'^-i - i In ( i^(^/2^)dl! ) for K<27r [RMT based interpolation - GOE ensemble] ^g) 

" 1 - iln((K/27r)''-i+l) for 2tt<k [RMT based interpolation - GOE ensemble] 

1 for 27r<CK [based on the diagonal approximation] 



The fact that we can deduce the explicit scaling functions (( p.27D and ( ^.28 )) for the classical correlations does not 



imply that their dynamical origin is understood. We propose some possible mechanisms in the next section. It should 
be also emphasized that the GUE scaling function is as little understood as the GOE scaling function, in spite of its 
apparent simplicity. See further considerations in section V. 

Let us assume for a moment that the scaling function possess actually a smooth crossover at K=2Tr. If this 
assumption were true then the tail of p{s) would be determined by the singularity at k ^ via the asymptotic 
relation C{k) ^ J^^p{s)ds that follows from ( 2.26| ). Therefore the asymptotic behavior of classical correlations 



would correspond to p{s) ^ 1/s'^. This argument concerning the asymptotic behavior oi p{s) is not-valid if a non- 
smooth crossover at k = 27r occurs, as suggested by RMT. Assuming that Kqm{k, L) actually obeys RMT prediction, 
one may perform an inverse Fourier transform of C(k) in order to find a complete expression for p{s) via ( ^.26|) . For 
d=2 assuming GUE statistics one obtains simply M] 

-^-9. N /sin(7rs)\^ ,„ „„x 

= . (2.29) 



(tts) 

This subject is discussed further in subsection III-F 



G. non-universal statistics 



For a generic ballistic billiard, due to level repulsion, the shortest quantal correlation scale is simply the average 
spacing Aq. It has been found that the short-range correlations on the scale Aq are well described by RMT. These 
quantal correlations are both generic and universal. If the billiard is characterized by a non-trivial structure, then 
one may find a shorter quantal correlation scale A_ due to splitting effect. (For example, one may consider the 
energy splitting in case of a double- well system). Such quantal correlation scale is neither generic nor universal. From 
semiclassical considerations one deduces that there should be also generic but non-universal quantal correlations. The 
latter should be found on large energy scales and are semiclassically related to the shortest POs. In order to deduce 
these long-range quantal correlations we should consider the basic statistical relation (2.14) in the CD regime where 
it reduces to Kqm{k,L) — KjjIL). The averaging window SL should be small enough in order to resolve individual 
POs, which is equivalent to the obvious requirement that i?gm(e) should not be multiplied by a cutoff window which 
is narrower than the corresponding correlation scale. The various correlation scales are illustrated in Fig.^. 
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(clustering) X 



FIG. 3. Illustration of the various correlation scales of the quantal spectrum (upper right axis) and of the classical spectrum (lower 
right axis). The relation of the generic non- universal correlations to either short POs or low-laying levels is schematically indicated. 

Similar considerations apply upon analyzing the classical spectrum. The actual length-spacing Ai ^ Lexp^—aL) 
has no statistical significance. The universal correlation scale Aq is much larger. Still, there may be a shorter classical 
correlation scale A_ due to e.g. clustering of POs. Such classical correlation scale is neither generic nor universal. A 
specific example will be discussed in section VI. There should be also generic but non-universal classical correlations 
on scales larger than Ao, which are related to the low- lying eigen-energies. These are deduced via the semiclassical 
relation Kci{k,L) = KQ{k) that holds in the QD regime. One should use a sufficiently small 5k in order to resolve 
these long-range classical correlations. This is equivalent to saying that the SCTF is capable of resolving individual 
energies levels of the quantal spectrum. 



H. Duality Within the Scattering Approach 



The classical dynamics of billiards can be conveniently described in terms of the Poincare maps M„ induced on 
appropriately chosen sections (ct) in the billiard phase space. A complete correspondence between the flow and the 
discrete dynamics is achieved when the section consists of the entire boundary surface and the tangential momenta 
at impact. Other sections, which are obtained by hyper-planes which cut through the billiard volume, and the 
components of the momentum in the hyper-planes may capture only a part of the billiard dynamics. 

The quantum mechanical analogue of the billiard Poincare map is the scattering matrix S{k\a) |^,^, which 
depends in a non trivial way on the particular section a and on the wavenumber k. In the semiclassical approximation 
the S{k) matrix is of finite dimension, which can be related to the volume of the d—1 dimensional section by 



Na{k) = 



kia 
TT 



(2.30) 



Where [x] stands for the integer part of x. It was realized long time ago that if the map induces chaotic dynamics 
on the section, then the statistical properties of the S{k; a) matrix are well reproduced by the predictions of RMT 
for the Circular Ensembles 45 1. In particular, 



1 

N 



(Itr^'f) 



1 

Ak 



fc+Afc/2 



\trS''{k;a)\^dk^Kf 



fc-Afc/2 



(2.31) 



Here, the integration is carried out over a k interval over which A^cr(fc) is constant N, so that k — 
7r((iV + l/2)i/('^-iVC) and Afc 7r/(2(d-l))iV-(''-2)/(d-i)/^^^ rpj^^ function kp{n/N) is known explicitly for 
the Circular ensembles (/3 = 1,2,4). See e.g. |4^. In ( |2.31 ) the k average replaces the ensemble average of RMT. 
This is justified since the correlator (iv[S{k)S^ {k-\-5k)\) vanishes once 5k exceeds the correlation length 5^ which is 
inversely proportional to the mean length of the classical orbit between successive encounters with the section a. For 
ergodic biUiard A/c » 5^. 

The semiclassical expression for trS'"(/c;cr) is 

tr^"(fc;a) « ^ 5f>"p(Ap)^ 



(2.32) 



where Pa-{n) is the set of primitive POs of Mg- of period rip which divides n so that 



POs which are 



conju gate by an exact symmetry are represented only once in Pa{n) and their multiplicity gp appears explicitly in 
( ^.32 ). The stability amplitud es (in clud ing t he Maslov and boundary indices) and the lengths are denoted by Ap and 
Lp respectively. Substituting ( 2.32 ) in ( 2.31 ) and separating the diagonal and non diagonal contributions we get 



n 



, n 



dx cos 



Ng + l/2' 



TTx Pa{x; n) 



(2.33) 
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where (g) is the mean multiphcity of symmetry conjugate orbits and Pa{x] n) is the classical POs correlation function 



Lp — Lpi \ sin[Afc(ip — Lpi)/2] 
Xa{n) ) [Ak{Lp-Lp,)/2] 



(2.34) 



where Xcrin) — Iryj-n}-^'^'^ is the classical correlation scale. The last factor limits the summation to pairs of POs 
whose lengths differ by (d-l)^^7V(''-2)/(''-i) at most. Equations ( |2.33H2.34| ) express the statistical duality between 
the spectra of eigenphases of the quantum S matrix and the classical POs of the corresponding map. Note that the 
correlation scale Xa{n) can be interpreted as the mean linear distance between successive impacts of the orbit on the 
section a. Therefore it is independent of a. It should be emphasize that the same result for A is obtained as in the 
analogous flow approach in spite of the fact that the POs are restricted to have a specific 'topological' period n. 

The construction above can be now used to further reveal the statistical properties of the classical spectrum. Let us 
denote by a the section which is built on the complete boundary, and by p any other section. Consider the POs of 
with period n^r- Only a fraction of them traverse p, and they belong to the set of POs of Mp with periods Up < Ua-. 
On the average na/ii~^ « '^p/^p~^ implying again that the classical correlation scale Ao-(n) is independent of a. On 
the quantum mechanical level, the scattering matrix S(k; p) will be of smaller dimension Np — Ncr{ip/£aY^^ ■ Thus 



(2.35) 



Let us compare the duality relations for the two sections. If and Up satisfy ( ^.35 ) then the left (quantum) hand 
side of (2.33) remains unchanged. The classical functions Pa{x] n„) and Pp{x; Hp) express the same statistical reality, 
but the latter uses only a fraction of the POs which contribute to the former. That is why pp{x;np) is rescaled by 
a factor (n^/up) which follows from (2.34). Thus, we can construct finer classes of POs which traverse smaller (but 
nested) sections, and which manifest the same POs correlations as contained in the original set of POs. However, this 
refinement is quite limited by the purely classical condition Xa-{n) <C i. The argument for this condition is as follows. 
Semiclassical considerations are valid as long as 27r <^ ki, else the scattering matrix becomes of dimension one. On 
the other hand, in order to deduce non- vanishing p(x; n) from the duality relation ( p. 33 ) the left hand side should be 
non-trivial, meaning N < n ot equivalently kX <^ 27r. Combining the two restrictions one obtains the fc-independent 
condition X^{n) <C i. Recalling the exponential proliferation of POs as a function of n, this condition implies that an 
exponentially large number of POs will be included in those classes. 



III. THE ORIGINS OF PO-CORRELATIONS: OBSERVATIONS AND SPECULATIONS 

It should be emphasized from the outset that we do not have a direct classical derivation for any of the results 
obtained above concerning the statistical two-point correlations of the classical spectrum. In this section we try to 
acquire, at least on a heuristic level, some insight concerning the origins of classical correlations. 

The first subsection introduces a simple statistical model for classical correlations using the notion of families. This 
family-structure is introduced on a purely formal grounds, without considering the question whether it is realized in 
actual dynamical systems. Still it provides an explanation for one of the puzzling features of the classical spectrum, 
namely, that its correlations are not apparent on the (exponentially small) level spacing, but rather on a much larger 
length intervals. 

In order to reveal the actual (statistical) structure of the classical spectrum we use various arguments. Some of 
the arguments constitute a variation of the quantal-classical duality argument for having classical correlations. In 
subsection B we introduce the notion of classification pl| ] . In subsection C we argue that one can find classifications 
that reflect the actual structure of the classical phase space. Thus it becomes plausible that some variation of a 
family-structure is indeed realized. Actually, we believe that a more complicated picture of hierarchical structure is 
appropriate (subsection E). The last subsection clarifies the role of resurgence in the theory of classical correlations. 

When considering POs of length Lj ~ L, one may define the geometrical scale Aq as the average distance between 
points where a typical PO of length L intersects a given surface of section. This scale appears naturally in the discussion 
of subsection C. There it appears as a limiting resolution scale for the classification of POs. In subsection D, we point 
that Aq also has the meaning of variation scale whenever parametric deformation of the billiard is considered. Then 
we argue that it should show up as a correlation scale in the two point statistics of the classical spectrum. Thus we 
obtain some heuristic understanding for the origin of classical correlations. 

We also try to speculate a formal approach for a future derivation of classical correlations. We suggest that an 
appropriate coarse-graining procedure may be used in order to approximate any generic billiard by a "graph" . This 
way metric correlations are transformed into a comparatively simple degeneracy structure. 
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A. families 



In order to gain some insight into the classical correlations, let us ignore the fluctuations of the instability amplitudes, 
and let us further assume that correlations between POs that correspond to different Maslov index are not important. 
These assumptions reflect our belief that the statistical correlations of the weighted classical spectrum are actually 
present also in the bare spectrum of lengths. To support this point of view one should notice that our result ( 2.22| ) for 



the classical correlation length Aq possess a very simple geometric interpretation: It is the average distance between 
neighboring points where a PO of length L bounced from the (d—l) dimensional billiard surface. 
Using loose notations, the two point statistics of the classical bare spectrum is 

S{x)—p{x;L) = ( ^ 5 (a; — (ij— Lpo)) ) —const (3.1) 

\ 3 / PO 

where the angle brackets denote averaging over reference POs whose length is Lpo ^ L, and const refers to the 
smoothed density. Note that the delta-function arises due to the diagonal term j — po (we assume no degeneracies) . 
If the rest of the lengths Lj are independent with respect to the reference PO, then p{x) ^ 0, which corresponds to 
A oo. However, if there is a small repulsion, then p{x) can be interpreted as the probability for finding a vacancy 
in the vicinity of the reference PO. 

A more detailed interpretation which seems to be relevant is as follows. Let us assume that the POs can be grouped 
into families such that each family constitutes a rigid spectrum ('ladder') by itself, with average level spacing A and 
strong repulsion. For clarity let us further assume that there are no cross-correlations between different families. If 
each family is characterized by the same two-point correlation function p{x), then also their union will have the same 
p{x). The latter interpretation makes it clear that the spectrum will look like Poisson on scales which are smaller than 
A. It is important to notice that the total density of POs grows exponentially with L due to their proliferation. Still, 
the average density of those POs that actually contribute to p{x) isp(O) ~ L^/'-'^^^\ It Implies that the overwhelming 
majority of periodic orbits j is uncorrelated with the reference PO. The rigidity of the spectrum manifests itself only 
on scales which are larger than A. 

This interpretation of classical correlations settles the apparent contradiction between the observation of Poisson 
statistics on small scales and the existence of non-trivial two-point statistics on the other hand. At this stage 
it is useful to draw an analogy with the statistical properties of a quantal spectrum that corresponds to a system 
with localization. The spectrum looks like Poisson on small energy scales, since neighboring levels correspond to 
well separated eigenstates. However, on large energy scales the rigidity of the spectrum is manifest. This rigidity 
reflects the fact that the spectrum constitutes a union of many local sub-spectra. Eigenstates that dwell in the same 
localization- volume constitute a family. 



B. classes 



In an attempt to understand the dynamical origins of the PO-correlations, it is useful to classify the POs according 
to some dynamical criteria, and to study the statistical properties of POs which share some common feature. One 
can consider for example a class of POs having the same number of bounces on the boundaries. A more refined 
classification can be obtained when symbolic codes are available. In this case one can classify the POs according to 
certain requirements on the codes. The most simple classification is obtained by gathering all the POs which are 
conjugate to each other by a symmetry operation such as time reversal or refiection. Such classifications may reflect 
certain statistical features of the classical spectrum. 

The classical spectrum is rigid since p{x) has a finite range A, and it satisfies the normalization J p{x)dx — 1. If 
the classification is arbitrary, then one will find that POs of a given class constitute a non-rigid spectrum. From 
now on we shall use the notion of 'class' and 'classification' in a more restricted fashion. By definition, POs of a 
class should constitute a rigid spectrum by themselves. A classification of POs into (mutually-exclusive) classes is an 
important tool for revealing the statistical properties of the classical spectrum. A classification is particularly useful 
if there are no cross correlations between different class-spectra. In some cases we can argue an a-priori existence of a 
classification scheme. At some other cases, as in section the classification can be guessed and verified a-posteriori. 
The existence of classes, with some restriction, has been argued in subsection II-H using the scattering approach point 
of view. In subsection C we shall argue that the notion of classification is appreciable to any billiard system. 

The most obvious example for useful classification of POs occurs while analyzing a disordered chain of cells. Such 
a complex billiard constitutes a quasi one dimensional diffusive system. POs whose period t is less than the ergodic 
time may be classified by the volume which they explore. POs that belong to different volumes cannot be correlated, 
since if the shape of one of the cells is modified, it will affect only those POs that explore it. 
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If the statistical model of subsection A were realized in actual dynamical systems, then the finest classification 
of POs would be into families. Further classification would not be possible. A direct understanding of classical 
correlations implies in some cases that a classification of POs into families has been indeed identified. An analytical 
theory requires further derivation of the two points statistics that characterizes these families, including the study of 
cross-correlations between families. It should be emphasized that different classifications may be found appropriate 
while analyzing non-universal rather than universal correlations. 

Assuming that cross-correlations between classes are negligible, it follows that the form factor can be decomposed 
in the following way 

K{k,L) = ^/^„(fc,L) = , (3.2) 

a a 

where Ka{L) is the diagonal sum that corresponds to POs of the a-class, while Ca{k,L) are the corresponding 
correlation factors. This formula will constitute the basis for our later study of complex systems. In case of a complex 
system, different classes may have different statistical properties. 



C. similarity 

If a the statistical model of subsection A were indeed realized in actual dynamical systems, then we would expect 
that POs of the same family would be geometrically related in some intimate way. In the present section we wish to 
use a variant of the quantal-duality argument in order to develop a notion of similarity. We do it by extending an 
idea due to Fishman and Keating pO[ . 

Consider a simple ballistic planar billiard. A magnetic flux line (j)Q penetrates the plane at a point vq which is 
located inside the billiard. The form factor K{k, L) with 0o 7^ will correspond to GUE statistics. Now, assume 
that additional flux lines (j)2, ■■■) penetrate the plane at the points (ri,r2,...). The form factor K^^y(k, L) still 
corresponds to GUE statistics, as well as its average K{k,L), where the average is over ((/>i,(/)2, ...). However, it is 
a straightforward exercise to prove that the semiclassical expression for the averaged form factor is given by 
where a denotes a class of POs that share the same set of winding numbers {vi, 1^2, ••■) with respect to the flux lines. 
Assuming that all the classes are char acterized by the same statistics. It follows that all the Ca{k,L) are identical 



with the GUE correlation factor ( 2.27 ). Thus, each of the classes constitute a rigid spectrum by itself. 

One may add more and more flux lines, thus achieving finer and finer classification of POs. Let us denote the 
typical distance between flux lines by i. Semiclassical considerations are valid as long as 27r ^ ki, where 27r/fc is the 
De-Broglie wavelength. Manifestation of rigidity is expected only in the quantal (Q) regime where fcAo(i) ^ 2tt. It 
follows that the semiclassical argumentation for rigidity of the a classes is restricted by the fc-independent condition 
Xo{L) ^ £. Thus we deduce that Xo{L), (which is essentially the mean distance between chords of a typical PO whose 
length is L), actually determines the similarity-resolution. 

It is plausible that in some sense, by using £ of order Aq, one will obtain the desired classification of POs into 
families. Literally, one may object this point of view using the argument that the locations (r'i,r2, ...) are arbitrary, 
therefore the classification into families is ill-defined. This objection can be dismissed by noting the analogy with 
localization theory (see end of subsection A). There, eigenstates may be classified by partitioning the physical space 
into 'blocks'. The partitioning into blocks is arbitrary. The limit of fine partitioning into small blocks, whose linear 
dimension is of the order of localization length, is problematic. Still, this limit is conceptually meaningful and leads 
to the correct identification of the family-structure. 



D. Bifurcations 



In this subsection we point out that the geometrical scale Aq has also the meaning of a variation scale whenever 
parametric deformation of a billiard is considered. We propose that this variation scale should show up as a correlation 
scale in the two-point statistics of the classical spectrum. We use quantal-classical duality in order to further support 
this conjecture. 

Consider the effect of some parametric deformation of a billiard system on the energy levels. We focus our attention 
on levels that are contained is some window around k. Changing a parameter b is considered to be a small perturbation 
if kn{b)—kn{b—0) <C Aq. Thus we can define the perturbative regime as 6 < be- For b ^ be the difference Sk = 
kn{b)—kn{b—0) is of the order Aq, and one expects to have gone via an avoided crossing. Thus we deduce that the 
correlation scale Aq can be identified with the variation scale 5k. The latter can be computed by extrapolating the 
leading order perturbative calculation up to the point where we expect that it will loose its validity. 
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Now we wish to speculate an analogous property for the classical spectrum. Consider POs that are contained in 
some window around L. The conjecture is that the correlation scale Aq should be identified with the variation scale 
6L. By definition, in order to determine the variation scale SL we should consider some parametric deformation of the 
billiard. Rather than having 'avoided crossings' we shall have 'avoided bifurcations'. The variation scale is calculated 
by applying linear analysis up to point where it is expected to loose its validity. Rather than considering a general 
deformation it is most convenient to analyze a specific deformation that takes us right to the bifurcation. Note the 
analogy between 'bifurcations' and 'crossings'. 

We consider a local deformation of the boundary. The normal distance of this deformation is b, and it involves a 
small surface area. We further assume that the linear dimension of the deformed surface is of the order Xo{L). At 
this stage Xo{L) is defined as the average distance between points where a typical PO of length L hits the surface. It 
follows that each PO of length L hits the deformed surface only once (on the average). If the deformation is small, 
meaning no bifurcation, then a PO keeps its identity, but its length is changed by Sj — 2bcos{6j) where 9j is the 
incidence angle (see appendix A for proof). The PO will disappear if b becomes of the order Ao(L). The deformed 
surface hits then one of the segments of the PO. Thus we conclude that the variation scale is Xo{L). It follows from 
our conjecture that this variation scale will show up in the two point statistics. 

It is interesting to analyze the significance of bifurcations from the quantal-classical duality point of view. Here 
we present a consistency argument that further substantiate the identification of the classical correlation scale A with 
the geometrical scale Xq. Considering the specific deformation of the previous paragraph, it is expected that the 5j 
are uncorrelated. The distribution of y — Si—Sj will be denoted by G{y). The average y is zero, and the standard- 
deviation is of order b. In appendix A we prove that if a classical spectrum is modified as in the present case, namely 
Lj Lj + 6j, then the new correlation factor is C{k) ^ [1 — G(k) ■ {1—C{k))]. Here G{k) is the Fourier transform 
of G{y). However, assuming that no symmetry breaking is involved, the new correlation factor should correspond to 
the same universal statistics as the old one. This will be true if either C{k) — 1 and/or G{k) — 1. The function G{k) 
has a crossover at k ^ 1/6, while the function C(fc) has a crossover at A; ~ 1/A, where A is the classical correlation 
scale (yet to be determined). Hence we obtain the fc-independcnt condition 6 ^ A. If the latter condition is not 
satisfied then our considerations will lead to inconsistency. On the other hand the only essential assumption for our 
considerations to hold is having no bifurcations. Therefore, it is plausible that the condition 6 <C A actually coincides 
with the condition b Xq leading to the identification of the universal correlation scale A with the geometrical scale 
Ao. 



E. Quantum Graphs and Shadowing 

Quantized graphs (networks) display spectral statistics which are reproduced to a high level of accuracy by the 
predictions of RMT H]. At the same time, the spectral density can be expressed in terms of an exact trace formula 
which involves summation over POs (loops). With an appropriate definition for the instability amplitudes, this trace 
formula looks formally the same as the SCTF of section II. Still, the graphs are sufficiently simple to enable a clear 
understanding of the origins of PO-correlations as well as their implications on the quantal spectrum |Q. 

A graph is a strictly one dimensional system. It is composed of a finite number of bonds (wires) which are connected 
at junctions (vertices). Classical dynamics is mixing because at each junction a particle takes a Markovian choice of 
the next bond on which it moves. The multiple connectivity of the graph induces the mixing nature of the system. 
The graph dynamics is Bernoulli, and as such it is characterized by an exponential proliferation of POs, to which 
positive Lyapunov coefficients can be assigned according to the probability to remain on the graph. The smooth part 
of the quantal spectral density is given by Weyl's law for a c? = 1 system. Thus the Heisenberg time is Lh = 2^4, 
where 4 are the lengths of the bonds (assumed to be non-commensurate to avoid non generic degeneracies). The 
oscillatory part of the spectral density is a sum of contributions from POs. Each contribution is a complex number 
whose amplitude is the square root of the classical probability to remain on the graph, and its phase is the action 
kLp where Lp is the orbit length. A topological phase factor e'^'^p plays a similar role to the Maslov index, /ip counts 
the number of time the orbit back-scatters from a vertex, (that is, the number of sequences of the type "aba" which 
appear in the code of the PO). 

The classical length spectrum for the graph has a simple structure - it is obtained by taking the linear combinations 
Lp = X^b'^fe^hi which are consistent with the connectivity of the graph, and are natural numbers. This structure 
of the length spectrum is revealed only if Lp > Lh = 2 ^ because for such lengths the orbit must be reducible 
to combinations of lengths of shorter orbits. We see again that the Heisenberg length Lh plays a natural role in the 
discussion of the morphology of the PO-spectrum. This is in complete accordance with the previous discussion of the 
duality concept. 

The fact that the Heisenberg length is constant (fc-independent) implies that one can generate the spectral two 
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point function using an arbitrarily large spectral interval (Sk). This way, one can reach the domain where the function 
Kqm{L) is composed of arbitrarily sharp spikes which resolve completely the length spectrum for lengths which are 
both smaller and larger than Lw. The smoothed Kqm{L) was compared with the predictions of RMT and was found 
to be in excellent agreement [Q. The only correlations of POs which survive this treatments are correlations between 
orbits which have the same lengths. One may consider this to be a complicated (generalized) case of degeneracy. Note 
that the instability amplitudes that correspond to degenerated POs should not have the same phase. For L < Lh 
there are hardly any length degeneracies, and therefore K{L) « L for L < Lh- For L > Lh the degeneracies of lengths 
increase rapidly with L, and the correlations between the POs which bring about the correct L dependence of the form 
factor are the correlations between the back-scatter indices /ip. This result demonstrates one important features of 
the POs correlations, namely, that POs carry both metric information (action integrals or just lengths) and topological 
information (Maslov indices and degeneracies). The correlations which are responsible to the quantum-classical duality 
sometimes reside in the metric information exclusively (Sinai Billiards) and sometimes in the topological information 
as was discussed above for quantum graphs. 

One of the most important recent developments in our understanding of the connection between classical chaos and 
quantum statistics was recently put forward by Bogomolny and Keating (BK) psf . They showed how to derive the 
non diagonal part of the quantal form factor by using information about short POs with Lp < Lh exclusively, and 
moreover, they assumed that these short orbits are statistically independent. We would like to discuss now the BK 
construction and show how it is related to the correlations which are discussed here. It is convenient to discuss this 
issue in the context of quantized graphs, since this system is both simple and exact. 

BK define an approximate spectral counting function N{k,L*) where the SCTF is truncated to include POs of 
length L < L*, and L* is yet to be determined. It is assumed that N{k, L*) provides a reasonable first approximation 
to the counting function, and it is improved by identifying the points fc* where N{k, L*) assumes the values n + ^ for 
all integers n. The resulting spectral density is 

d{k) = Y,5(N{k,L*)-n-]\^-^^fl = d{k,L*) (_i)™e2""^('=.i*) (3.3) 

n ^ ^ m— — 00 

The main purpose of this construction is to ensure the 5 function structure of the spectral density. This is done, 
however, at the cost of producing a spectral density which might differ from the real one in detail. This can be easily 
understood when ( |3.3| ) is expanded and the length spectrum is deduced. The length spectrum beyond L* consists 
of all the sums 'Y^rripLp where here the vrip are integers, and the basic lengths Lp < L*. There are no restrictions 
on the combinations of the Lp, which for quantized graphs implies that the rules of the symbolic dynamics are not 
obeyed! In other words, even the topological entropy of the synthetic POs is not that of the original system. BK used 
their procedure to obtain averaged spectral correlations which might be less sensitive to the problems inherent to the 
method. 

Using the intuition we derive from the work on quantized graphs, we would like to propose a heuristic understanding 
of classical correlations in billiard systems with d > 1. For billiards Lh is fc-dependent as in Figjlj Suppose that we 
need to perform a classical summation of the type exp(ifcLj). We wish to approximate the true lengths Lj by 
some new set of lengths that corresponds to a coarse grained phase space. Given k one can define a spatial resolution 
scale 27r/fc. It is plausible that this coarse graining scale is appropriate in order to obtain a faithful approximation for 
the classical sum. However, we do not have a purely classical argumentation for this speculation. If this speculation 
is true, then we may approximate the length of long POs by a sum of the form X^b^-fc't" where {lb} is a set of 
comparatively short POs that form a skeleton structure. The idea that long POs are supported by short POs is called 
'shadowing' |^6|. We can estimate the sum J2b for the skeleton which is required for a faithful representation of all 
the long POs of the coarse-grained phase space. We do it as follows: due to the coarse graining each PO of the skeleton 
has a thickness {2n/kY~^. The total volume of the skeleton structure is X]h(27''/^)'^~^4, and it should be equal to Vl. 
The billiard is now approximated by a "graph" . The Heisenberg time for this graph is Lh = 2 l^ ^ flk"^"^ , which 
is the generic quantal result. 

Using the shadowing concept we can therefore obtain a heuristic understanding for the origin of classical correla- 
tions. By improving gradually the coarse-graining resolution of phase space, the hierarchical structure of the classical 
spectrum is revealed. POs that are identical in length on a given resolution scale, will become a rigid family if the 
classical spectrum is probed with better resolution. 

F. Resurgence 

It turns out that the statistics of composite periodic orbits (CPOs) is simpler that the statistics of primitive periodic 
orbits (PPO) [^. The CPO-spectrum -^(n.ra,...) ~ J2''^j^3 is obtained from the PPO spectrum, where the rj are 
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integer repetitions. There is a strong analogy PPO <S4> primes, and CPO <S4> natural-numbers (see section IV-F). 
Actually, it turns out that the two point statistical correlations of CPOs is essentially the same as that of the natural 
numbers, meaning that the scaling function is universal and d-independent. Namely, 

TTS 

The analogy to our formulation has been extended |2^. The role of the density Pqm{k) ^ A:„) is taken by 

the Zeta function ^ n(^~^n)- The two points correlations of the former are Rqm{k,e), while the two points 
correlations of the latter are C'qm{k,e). The latter auto-correlation function is essentially real and therefore the 
corresponding form factor Kq^{k,L) is symmetric, namely satisfies the functional relation 

K''{k,L) = K%k,LHik) - L) . (3.5) 

The quantal form factor Kq^{k, L) is related to the two points statistics R^i^'~'{x, L) of the CPO spectrum. Therefore 
( |3.5[ ) can be re-interpreted as a classical relation between the statistics of long POs and the statistics of short POs. 
The existence of such relation is hardly surprising, it should follow from the observation that long POs are supported 



by shorts ones. A classical derivation of (3.5) would be a major step in the understanding of classical correlations. In 
p2| the functional relation has been re-expressed in terms of PPO-statistics, leading to the well known relation [^3|j25[ 
between the diagonal behavior of K{t) for short times, and its 'off-diagonal' behavior in the vicinity of Heisenberg 
time. 



IV. APPLICATIONS, GENERALIZATIONS AND EXAMPLES OF DUALITY 

The formulation of quantal-classical duality will be extended to the case of either general chaotic system (subsection 
A) or general integrable system (subsection C). We also discuss specific adaptations of the general formulation: The 
introduction of magnetic field and the analysis of periodic chains (subsection B); The theory of the 3D torus billiard 
(subsection D); The theory of the 3D Sinai billiard with mixed boundary conditions (subsection E). Finally (subsection 
F) we discuss, using the same approach, the statistical properties of the primes. 

Both, the classical spectrum of the cubic 3D torus billiard, and the 'classical' spectrum of the primes, are char- 
acterized by a non-trivial correlation scale. These two examples, together with the generic case of chaotic billiard, 
demonstrate our claim that there is no a-priori relation between the classical correlation length A and the average (or 
typical) spacing between actions AL. Indeed, if we summarize these three cases we get: 

A <C AL for the primes 

A ~ AL for the cubic 3D torus 

A ^ AL for a generic billiards . 

Note that both A and AL are functions of L. These functions have no a-priori inter-relation. 



A. General Chaotic Systems 



For general time-independent Hamiltonian system the SCTF is: 



Xn 



2TTh 



t r ^\det{M;^I)\ 



exp 



rSp{E) 



(4.1) 



where {En\ are the eigenenergies. The actions Sj = rSp correspond to the primitive POs and their repetitions (both 
positive and negative). The primitive periods are Tp, and the Monodromy matrix is Mp. The effective Maslov index 
Uj — rvp is incorporated with the appropriate (negative) sign that corresponds to time-domain FT conventions. 

In order to apply the formulation of section II we may use one of two strategies. The first strategy is to consider 
the relation between the ?i-spectrum and the S'-spectrum. The former is obtained by fixing the energy E and varying 
h. It is convenient to define the quantal variable a = l/?i, and to cast (4.1) into the following form 



— 6{a—an) = Aj exp (iaSj) 



(4.2) 
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where w„ = a{dEn/da) are the so-called velocities. The dependence of the actions Sj on the energy E is not emphasized 
in the above formula since the value of E is fixed to some constant value. With this notations the formulation of the 
relation between the quantal a-spectrum and the classical S'-spectrum is exactly as for billiards, the only difference 
being that now also the quantal spectrum is weighted. 

If one insist to fix a and to consider the statistical properties of the _B-spectrum, then a different strategy is 
required. Picking an energy E, one should assume that it belongs to an energy interval [Ei,E2] where the POs are 
structurally stable (See remark |^), meaning that there are no bifurcations of POs while E is varied. Thus Sj{E) is a 
well defined function in this range. The time periods are Tj — {dSj{E)/dE). The general statistical theory of section 
II applies with small modifications. It is also essential to make the approximation Sp{E) = Sp{E) + {E—E)Tp, and to 
define a scaled energy variable co = {E—E)/h. Note that {E2—E1) is classically a very small energy interval, while 
ujmax = {E2—Ei)/h IS scmiclassically very large. The SCTF is now 

y 27r(5(aj— tj„) = N exp (iaS'j + iwTj) , (4-3) 

^ — ^ osc ^ — ^ 

j 

The Two point correlation function Rci{a,e) of the classical spectrum will be considered, where e corresponds to 
difference in the energy variable uj, and w = without loss of generality. In order to establish the basic statistical 
relation Kqm{a,t) = Kci{a,t), it turns out that the proper definition of the classical two point correlation function 
should be 

Rci{x.t) = (^A,Aj 5{x~{S,~S,)) 5{t-Ti)^ - const , (4.4) 

where const is the subtracted smooth component. This two point correlation function docs not correspond to any 
one-point density unless there is a scaling relation such that Sj oc Tj. Actually it is sufhcient to assume a strong 
statistical correlation within families, rather than an exact global scaling relation. 

There is a better formulation of the general quantal-classical duality, which enables to treat on equal footing both 
the classical and the quantal spectra, as in the special case of billiard systems. This formulation is based on a time- 
domain version of the SCTF. The validity of standard FT procedure is somewhat more subtle, and therefore we shall 
only sketch the main idea. The time domain version of the SCTF is 

1/2 

. En \ 1 / dE„ \ t ( . rA-D . TT 



where Aj = rA-p denotes the ]uU action (integral of the Lagrangian) of POs whose period is t. With the appropriate 
definitions, this relation can be cast into the form 



27r 



1/2 



E 



exp(— itct;„(Q:)) = Aj exp(m^j(t)) , (4.6) 



where a = Xjh. An optional way to write the same is 

^T{i^'^t\Pqm{(^,'^) = ^T^{A'^a\Pcl{A,t) . (4.7) 

Note that the classical density pdiA^t) depends parametrically on t, while the quantal density pqrn{ct,uj), depends 
parametrically on a. Note also that the smooth a dependent factor (277/0;)^/^ of equation ( |4.6| ) should be absorbed 
into the definition of the quantal density. Both sides of the equation are functions of (a,i). Squaring them, and 
performing statistical averaging, one obtains the basic statistical relation Kqm{o:,t) = Kci{a,t) in complete analogy 
to the billiard formulation of the quantal-classical duality. 



B. Billiards with magnetic field, Periodic chains 



For a billiard in a miiform magnetic field one may use the general trace formula (4.3) with a replaced by k, and Tj 
replaced by Lj, and Sj — Lj + bAj. Here Aj is the net area that is enclosed by the PO, while b — (qB / c) / (mv) is the 
scaled magnetic field. For simplicity units of length are used also for the action. If the magnetic field is concentrated 
in one flux line then 

5, =i,+H (4.8) 
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where Vj is the winding number and h = {q4> / c) / (mv) is the scaled magnetic flux (Note that it has now the dimension 
of length). In particular one may consider Aharonov-Bohm ring geometry. The same formulation applies. 

Periodic chains can be treated using the same formulation that is applied for the semiclassical study of Aharonov- 
Bohm ring. The quasi-momentum q is a constant of the motion, and therefore can be regarded as a parameter. 
Eigenstates that correspond to q should satisfy the Bloch condition 'ip{x + a) = eKp{iqa/h)ip{x). This is equivalent to 
solving one-cell problem with periodic boundary conditions and a scaled magnetic flux b — {qa)/{'mv). The winding 
number i' of the folded PO characterizes the periodicity of the corresponding unfolded PO. 

The form factor K{b,k,L) depends paranietrically on b. If one averages over &, then one obtains K(k,L). It is a 
simple exercise to prove that the corresponding classical two-point function satisfies 

K,i{k,L) = ^i^,(fc,L) = J2c,{k,L)K,{L) , (4.9) 
meaning that cross-correlations between classes of POs that are distinguished by their winding number should be 



ignored in computation of b averaged two-point statistics. The last equation would become a special case of (3.2) if 
cross-correlations were negligible. 



C. General Integrable Systems 

For general time-independent integrable Hamiltonian system with d degrees of freedom the SCTF is: 

E^(^-^«) ... - ir;kmY.Y.^v^r^''^'-^A"-^~^ (4-10) 



Here the sum is over rational tori p, whose frequency vector is = (27: /Tp)rrip, where nip is an integer vector of 
relatively prime integers and Tp is the primitive period. The scalar curvature of the torus is Kp, and the Maslov-Morse 
indices are a and (3. Using the same strategy as in subsection A, and using analogous notations, we can cast this 
formula into the form 

(d^i)72 ... " ^Ajeyi^{iaSj + iujTj) , (4.11) 



where a = 1/h. The quantal form factor Kqrn{ct,t) that correspond to the left hand side, should be equal to the 
classical form factor Kci{a,t) that corresponds to the right hand side, defined as in subsection A. 

The asymptotic behavior of the quantal form factor in the deep quantal domain (Q) is determined by the quantal 
diagonal approximation, namely Kqm{t) « Kq^o). It reflects the self-correlations of the levels, and therefore, if we 
assume no degeneracies, it is determined merely by Weyl law. One obtains 

Kq,n{a,t) « KQ{a) = \ ^^j^^ deep quantal regime. (4.12) 

[Zn ) ah/ 

Similarly, in the deep classical regime (C) one may apply the classical diagonal approximation. Here one should utilize 
the version of the Hannay and Ozorio de- Almeida sum rule that applies to integrable systems . See also appendix 



B of |27|. Assuming no degeneracies one obtains 



Kci{a,t) w Koit) = /o L 1 ^^}^^ deep classical regime. (4.13) 

(2n)"-^^ ah 

Quite surprisingly KQ{a) does not depend on a, Koit) does not depend on t, and both are equal. The quantal 
diagonal approximation is valid in the deep Q-regime and the classical diagonal approximation is valid in the deep 
C-regime of the (a, t) plane. It is therefore very tempting to assume that the form factor K(a, t) is constant over the 
entire (a, t) plane, and it does not experience any actual crossover. We shall show below that this is indeed the typical 
case: For a generic non-degenerated spectrum that corresponds to a generic integrable system both the quantal and 
the classical spectra are of Poisson type. 

The time domain formulation of the quantal-classical duality can be extended also to the case of integrable systems. 
Rather than considering actions Sj of rational tori on the energy surface E, one should consider the actions Aj of 



rational tori with a period t. Then one obtains a SCTF that can be cast into the form of (4.6) with the h dependent 
factor (27r/a)^/^ replaced by (2-k j aY"!"^ ^ where d is the number of degrees of freedom. 

The Poissonian nature of both the quantal spectrum and the (dual) Poissonian nature of the classical spectrum, 
can be deduced using the following argumentation. The number of energies up to an energy E equals essentially to 
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the number of lattice points / = hn, such that E{n) < E. Analogously, the number of actions is given by the number 
of lattice points uj = (27r/i)m, such that A{7n) < A. The quantal and the classical counting functions grow by 1 
whenever the respective constant E surface or the constant A surface cross a discrete lattice point. The crossings 
occur at values of E (or A) which depend in a very complex way on the properties of the moving surface, and introduce 
the spectral fluctuations which decorate the smooth counting functions. Berry and Tabor where the first to relate 
the spectral statistics of integrable systems to the problem of counting the number of lattice points within a surface. 
They have shown that the fluctuations for generic surfaces are Poissonian. From this point of view, the constant E 
or constant A surfaces are similar, and therefore, the statistics for both counting-functions must be Poissonian. The 
argument should be modified somewhat if one considers the actions Sj rather than the actions Aj . Each point in the 
space labels a torus, and the energy surface is represented by some closed surface. Tori that correspond to the same 
frequency ratio m are located along a ray. The intersection of the ray with the energy surface determines the period 
and the corresponding torus is such that uj = (27r/i)m. Let us associate with the ray a length that correspond to the 
action Sj of that rational torus. Locally the tips of the rays will form a rectangular grid structure in any direction. 
However, globally they form an irregular array of points. The number of actions Sj < is equal to the number of 
tips that are contained in a sphere whose radius is S. Rather then having an irregular surface that moves in a regular 
grid we have a 'regular' sphere that moves in an irregular grid. Therefore the same considerations that lead to the 
deduction of Poisson statistics apply also in this case. 



D. Example: The 3D Torus Billiard 

In this subsection we discuss the 3-dimensional torus billiard. It will serve as an interesting illustration of the duality 
relations presented in the previous subsection. We discuss first the general 3D torus with unequal edges (01,02, as). 
The quantal spectral density is given explicitly as: 

dq^{k) =Y,Hk- M , % = 27rv/(pi/ai)2 + (p2/a2)' + (ps/aa)' ■ (4.14) 
The classical density is: 

dciiL)= ^J(L-Lp-), Lp- = V(Piai)' + (-0202)2 + (psas)^ • (4.15) 

The Berry- Tabor semiclassical trace formula for the quantal density of states is obtained by performing Poisson 
summation of dqm which can be carried out exactly: 

, , VP Vk sm(kLs) -, , Vk f+°° , , sin(fcL) 

In order to cast ( |4.16|) into the explicitly dual form pqj^ = TT pd one should define 

2-K iV 
Pqm{k) = — K„(fc) + dq,^{-k)]^^^ , p,i{L) = — [d,i{L) + dci{-L)] . (4.17) 

From these expressions we observe that both the quantum and classical spectra are supported by 3D rectangular 
lattices, and hence they should have the same statistics. The argument at the end of the last subsection further suggests 
that this statistics is Poissonian. Moreover, both the classical and the quantal asymptotes (diagonal approximations) 
of the form factor have the same value: 

KD{L)=KQ{k) = -, (4.18) 

TT 



which is indeed independent of both k and L, confirming (4.12, 4.13). All these arguments strongly suggest that the 
form factor is a constant for all k,L: 

K{k,L) = —. (4.19) 

TT 

We note that the conventional form factor which is calculated from dqm rather than from is Vk'^/ir and is not 
purely classical for small L values (classical diagonal regime). 

We now turn to the special and interesting case of the cubic 3-torus: 01=02=03 = o. In this case, both 
the quantum and the classical spectra are directly proportional to the vectors of the cubic lattice p <E '2? . The 
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vectors p = {l,m,n) exhibit a geometrical degeneracy of their lengths p = + w? + n? due to permutations of 
the indices and/or sign changes. Much more interesting is the number-theoretical degeneracy of 2^, for example 
5^ + 6^ + 7^ = 1^ + 3^ + 10^ = 110. Denoting by g{p) the overall degeneracy of p: 

g[p) = i^{K^7l? : \R\^p] , (4.20) 

one have the following results : 

^^^=/3p, = const. (4.21) 

{gyp)) 

Applying this result to the classical and quantal diagonal approximations give: 

KoiL) = ^L, KQ[k) = -^k. (4.22) 

These asymptotes are in sharp contrast with the Poisson (generic integrable) case and are actually reminiscent of RMT 
predictions (generic chaotic) in 2 dimensions! This unexpected behavior of "chaotic" form factor for an integrable 
system is due to the number-theoretical peculiarities of the cubic 3-torus. This is an interesting example, when the 
classical-quantal duality is explicit, and yet a statistics which is more interesting then Poisson is obtained. 
The expected break time L* occurs when the two asymptotes intersect: 

L*{k) = ^k. (4.23) 
This equation defines the crossover line in the (fc, L) plane between the classical and the quantal regimes as discussed 



in (2.21). Thus we deduce that the classical correlation scale is given by: 



\{L) = - (4.24) 
while the quantal correlation scale is: 
2tiV 1 



A(fc).^-J-. (4.25) 

We can relate the scale A to the minimal distance between lengths of POs. Indeed, due to the structure of the POs 
Lp- = ap = a\/intcgcr, the minimal difference between lengths near length L is 

AL^^ (4.26) 
and consequently 

A(L) = 2AL . (4.27) 

This indicates that the breakdown of the diagonal approximation is related to the "number-theoretical rigidity" of 
spacing 1 between integers. More over, it indicates that there is only one underlying family of POs in the sense that 



was discussed in subsection III A . Analogously, we get for the quantal scale: 

A(fc) = 2Afc , (4.28) 
where Afc is the minimal spacing between consequent energy levels. 

E. Mixed boundary conditions and Sinai billiards 

We present in this subsection a class of billiard systems which will be central to our numerical studies. The most 
general boundary conditions under which the quantum billiard problem is self-adjoint can be written as p^ : 

K cos l3 ip{r) + sin P dftipiT^) = , r G (boundary of the billiard) , n = normal pointing outside , (4.29) 
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where the angle /3 interpolates smoothly between Dirichlet 
with the dimension of a wavenumber. It was shown in ||3 
modified to incorporate mixed boundary conditions: 

N{k;K,f3) = ^e[A;-fc„(K,/3)] 

n 

= N{k; K, /?) + ^ Aj exp |j 



= 0) and Neumann (/3 = 7r/2) cases and k is a constant 
that the Gutzwiller trace formula for billiards can be 



PO 



kLp '2^P 



$p(k,/3) = ^0(6'f;K,/3) , 0=(-2)arctan 



KCot/3 
k cos 6 



(4.30) 
(4.31) 



In the above, A, 



l/{2nir\J\AIp^ — I\) are the instability amplitudes, p is the index of primitive POs, r e Z is the 
repetition number, h'p is the Maslov index due to focusing only. Hp is the number of reflections from the billiard's 
boundary of the p'th primitive PO, and 0f are the reflection angles with respect to the inwards normals. Explicit 
expressions for N{k; k, (i) appear in js^. The parameter n can be generalized to become a piecewise constant function 
on the boundary, with the obvious modifications of <I>p. This is a very powerful tool to selectively filter semiclassical 
contributions of POs. If we change the boundary conditions only on a segment of the billiard (e.g. the disc in the 2D 
Sinai billiard) then orbits are selected according to the number of reflections from this segment only. Let us proceed 
by taking the derivative of ( [I.3l[ ) at /? = (near Dirichlet) and get the expression 



^N{k-K,P) 



/3=0 



S{k - kn) 



/3=0 



d_ 

dp 



N{k;K,f3) 



/3=0 



PO 



(4.32) 



fc„ = fc„ (Dirichlet) 



2irk ' 



i=l 



cost 



Equation (4.32) is similar in structure to the standard trace formula, but for a few modifications. The quantal part is 
essentially the density of states for Dirichlet boundary conditions dressed by the "velocities" Vn = 9fc„(K, /3)/9/3|/3=o- 
The phases in the semiclassical summation over POs are the same as for the Dirichlet case, and the semiclassical 



amplitudes are modified relative to the standard case: Aj 



A,B,. 



The most important application of this formula 



is to efficiently filter the non-generic contributions of the "bouncing-balls" in the 2D and 3D Sinai billiards. The 
bouncing-balls are families of neutral POs which are of classical measure 0, whose effects completely shadow the 
contributions of generic chaotic POs in the semiclassical analysis p3|-^. Instead of subtracting their contributions 
explicitly which is a difficult task in 3D, wc use (4.32): Taking the derivative with respect to /3 is essentially a 
subtraction of two spectra with different boundary conditions on the sphere. Since the bouncing balls never reffect 
from the sphere, their leading contribution is trivially eliminated. Moreover, Bj oc ^iCOsOj = if and only if the 
PO is wholly tangential to the sphere, thus we eliminate also the closure s of t he bouncing ball manifolds (see figure 
which can otherwise overwhelm the generic contribution. Equation (4.32) therefore relates the spectra of Sinai 
billiards with contributions of generic periodic orbits only. 
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Finally, to cast (4.32) into the standard from Pqm{k) — TTpci{L) we divide both sides by k: 

^ u„(5(/c - A;„)|„^^ = ^ JLjexp(ifcLj) , {;„ = ^ , Aj = exp {-i^Vj - injir) . (4.33) 

po fc„ /c v ^ / 

The left hand side is a weighted density which is purely quantal, and it is a Fourier transform of a purely classical 
density J2j AjS{L — Lj). 

F. Theory of the primes 

It is instructive to test our formalism in the special limit "d 1" which corresponds to the spectrum of the 
primes Q. The "classical spectrum" consists of the real numbers Lj = rlog(p) where p are the primes and r is 
the repetition index. The "quantal spectrum" is defined by the non trivial zeros of the Riemann Zeta function, 
namely C(l/2 + i^n) = 0. By the Riemann hypothesis the spectrum {fc„} consists of real numbers. The "quantal" 
spectrum is related to the "classical" one by a trace formula which is formally identical with the SCTF ( |2.5| ) with 
Aj = Lp exp(— rLp/2). The average density of the classical spectrum is (^5{L—Lj)) — exp{L)/L and therefore the 
Haney-Ozorio sum rule (^) is satisfied. The average density of the quantal spectrum is {J2 2TT5( k—k n)) = log(A;/27r). 
This functional form corresponds in some sense to the dimensionality From here (via ( ^.21 )) it is expected 



that the classical spectrum is characterized by a correlation scale Aq = exp(— L). This correlation scale expressed in 
natural units (via the relation L = log(n)) is simply An = 1. Thus Aq has a trivial number theoretical interpretation. 
It simply reflects the fact that the classical spectrum is supported by an underlying grid that is formed by the natural 
numbers. The manifestation of the Aq scale in two-point correlations in this case is not a "conjecture" but rather a 
simple mathematical fact. 

Rather than going on "backwards" as in section II-C, let us use at this point the so called Hardy-Littlewood 
conjecture. Following Keating [Q the statement is that if n is a prime, than the probability that n+An is also a prime 
is reduced, the factor being (1— l/(2An)). The statement holds asymptotically, but we shall assume that it actually 
holds for any l<An. Thus the two point correlation function of the classical spectrum is characterized by repulsion 
with p{x) — (l/(2An)) • with x = Alog(n). It follows that that p{x) = l/{2Lx) with lower cutoff at 

A = exp(— L). 

Having "derived" an expression for the classical correlations, we should verify now that indeed the basic relation 
Kqrn = Kci IS Satisfied. To be more precise, we wish to examine whether using p{x) one can calculate semiclassically 
the departure from the diagonal approximation, and whether the correct asymptotic behavior Kqm{L) — )■ KQ{k) 
is achieved. Else, one should conclude that the trace formula would not converge to a density that corresponds to 
the spectrum of Riemann zeros^ Obviously, for k <C 27r/A the diagonal approximation C{k,L) « 1 is valid. For 



27r/A ^ k one easily find via (2.26) that C{k,L) « log(fc/27r)/L. Hence K{k,L) Ko jk) as is expected. Note 
that the normalization of p{x) is satisfied in an asymptotic sense, and is reflected via ( 2.26| ) by the observation that 



C(A;, L) in the deep quantal regime (large L) 

The spectrum of the primes is an example for system where the correlation scale A is smaller than the average 
level spacing. This should be contrasted with d-dimensional systems where the mean level spacing is very small due 
to the fact that the spectrum is apparently a union of many independent sub-spectra (families). This observation is 
important since it demonstrate that 27r ^ /cAL is neither sufficient nor necessary condition for the validity of the 
diagonal approximation | Aj exp(iA:Lj)p « In case of the primes A ~ exp(— L) which leads to L ~ log(fc). 

If one, by mistake, considered AL Lexp(— L) as the relevant scale, one would deduce a longer breaktime (still it 
will be roughly logarithmic in k). The same wrong time scale would be deduced in case of billiards. However, for 
billiards it would be much shorter than the actual Heisenberg breaktime, as discussed in section II-E. 

The non-trivial statistical properties of the classical spectrum, in case of the primes, are implied by the fact that the 
spectrum is supported by an underlying grid that is formed by the natural numbers. In analogy, one may speculate 
that the non-trivial statistical properties of the classical spectrum, in case of billiards, is due to the 'fact' that long 
POs are supported by short ones (nGl. 



V. COMPLEX BILLIARDS 



We use the term 'complex billiard' when the chaotic dynamics of the billiard is characterized by a non-trivial 
approach to ergodicity. Such 'non-universality' should have manifestations in the one-point and in the two-point sta- 
tistical properties of both the classical and the quantal spectrum. It is associated with having non-trivial classification 
of POs. As a first example we mention the 3D Sinai billiard. We shall see (section VI) that it is essential to classify 
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the POs by the number of times that they are reflected by the inscribed sphere. In the following subsections we 
discuss the complexity that arise due to either spatial proximity or global symmetry. 

It is important to understand that the existence of 'classical correlations' is not restricted to simple ballistic billiards. 
On the contrary, the effect of PO-correlations should be much more significant when complex systems are involved, 
since relatively short breaktime-scales may show up. Our interest in complex systems is twofold. On the one hand 
we want to develop a semiclassical theory which is not limited by the diagonal approximation. On the other hand we 
want to test the limits of the semiclassical approach. There is a long standing question concerning the feasibility of 
having PO-theory for the spectral properties of system with either band structure or Anderson localization. We want 
to use our present level of understanding of the SCTF and of classical PO-correlations in order to give an affirmative 
answer to this question. For clarity of presentation we shall sometimes use the symbol t rather than L in order to 
denote lengths. This is done in order to adhere to common notations and terminology. 



A. Symmetries and Symmetry Breaking 

If a system is characterized by a symmetry, for example a reflection symmetry, then it is natural to consider first 
the POs of the folded system, and then, from them, to generate the POs of the unfolded system. Most of the POs of 
the unfolded system will not be conjugate to themselves under symmetry. Therefore, they will have some degeneracy 
g. One may regard this degeneracy as the simplest example for a classification of POs. The corresponding two points 
statistics will be characterized by a correlation factor 

C{k,L)=gCo{kXg{L,n)) . (5.1) 

It is assumed here that further classification of POs is not required, else an explicit class index a should be introduced. 
The scaling function Co refer to the bare length spectrum which is obtained by ignoring the degeneracy. The scaling 
constant should correspond to the folded billiard, therefore Xg{L,fl) = Xq{L,VI/ g). It turns out that the latter 
scaling strategy is useful also in the case of time reversal symmetry, though the argumentation for having Q./2 is 
less transparent. As a general rule phase-space volume rather than physical volume should be relevant. The scaling 
functions Cq{h) for GUE statistics as well as for GOE statistics has been already displayed in section II-F. 

Next, we should discuss the effect of symmetry breaking due to some deformation of the billiard. We are still 
interested in the two-point statistics of POs whose periods are Lj ^ L. We assume a small deformation h of the 
surface, that does not involve any bifurcations of them. Accordingly, the POs change their length Lj Lj + 5j 
because of the deformation. Our main concern is the statistics of the differences y = Sj—Si. Due to the central 
limit theorem we expect Gaussian statistics with dispersion y/VL where V is some constant that formally resembles 
a diffusion coefficient. It is convenient to define the deformation parameter via the relation D — /Co, where Cq is 
the mean chord (linear dimension of the billiard). The parameter 6, which has the dimension of length, constitutes 
the natural measure for the surface deformation. With these notations the distribution of y is 



G{y)^^===^e^p{-,^^y^L__) . (5.2) 



1 ( y_' 

^2niL/Co)b ""^K 2iL/Co)b 
Besides the universal correlations of the Lj , one may define a non- universal correlation scale 



X*{L) = ^WUjb. (5.3) 

Unlike the universal correlation scale Ao (L) , the non- universal correlation scale A* (L) is a monotonic ascending function 
of L, as illustrated in Fig.^. Following the discussion in section III-D it is conjectured that the validity of (5^) is 



restricted by the condition A* (L) <C Ao (L) . This determines a classical length scale Lc and a corresponding time scale 
tc beyond which the effect of bifurcations should be taken into account. It is natural to associate this time scale with 
the crossover Ag — *■ Ao of the universal correlations. From now on we shall use the notation t rather than the notation 
L in order to denote lengths. Accordingly, we shall use the notion of 'time' as a synonym for the notion of 'length'. 
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FIG. 5. The universal correlation scales Aq and Xg and the non-universal correlation scale A* are plotted against t. The horizontal line 
facilitate the determination of the various timc-rcgimcs for a given k. 



Recall that the Heisenberg time tnik) is determined by the condition kXo{t) ^ 2tt, (See section II-C). litnik) ^ tc, 
then there will be no significant manifestation of the symmetry breaking. On the other hand, if tc <C tuik) as in Figj^, 
then there will arise a non- universal breaktime t*{k) that is determined by the condition kX*{t) ~ 2tt. It should be 
emphasized that unlike the universal breaktime scale (tn), the additional breaktime t* does not signify a crossover to 
a quasiperiodic behavior in the present context. (This should be contrasted with a situation which will be discussed 



in subsection E). The correlation factor (5.1) is modified due to the symmetry breaking as follows (see appendix A) 



C{k, L) = 1 + (.g-l)-G(fc) [Classical regime, where Co(fc, L) = 1] . (5.4) 
Above G{k) = exp{-t/t*) is the Fourier transform of G{y), and t*{k) = 2Co/{bkf 



We wish now to discuss whether there are circumstances in which a more complete version of (5.4) exists, which 
describes the actual crossover which is associate with the symmetry breaking. Taking into account the similarity of 
statistically related POs, one should specify whether the deformation preserves the integrity of families. If the answer 
is negative, then the 5j may be considered uncorrelated random variables, and then we have as in section III-D the 
result C(fc, L) = [1 + G{k)-{gCo{kXg) — 1)]. However, as discussed in section III-D the validity of the latter expression 



should be limited to the regime where it coincides with (5.4) . However, if integrity of families is preserved, then one 



may assume that Sj for POs of the same family are the same, but otherwise they are uncorrelated random variables. 
Then one obtains 

C{k,L) = [1 + {g-l)G{k)]-Co{kXg) (5.5) 

The validity of the latter formula is not restricted by bifurcations, since it has been assumed that the family structure 
is essentially preserved. Still, the global description of the symmetry breaking is not satisfactory since the Xg Xq 
crossover is missed. The actual reason for this failure is further illuminated in the next paragraph. 

Breaking of time reversal symmetry deserves a particular attention. It is simplest to discuss Aharonov-Bohm ring 
geometry, since changing the magnetic field does not involve bifurcations of POs. Rather, Lj —^ Lj +S_j, where 



Sj — bvj. Here the POs are classified by their winding number v and 6 is the scaled magnetic field (see (4^)). The 
distribution of v is known to be Gaussian with dispersion \/'DL/ a^j, where T) is the coefficient for diffusion along the 
ring (one should imagine an unfolded representation of the ring as a periodic multi-cell structure) , and ao i s th e length 
of the ring (length of unit cell in the unfolded structure). The distribution of y is accordingly given by (^^), where 
£o = Qq/-^/'^- Each winding- number class constitute a rigid spectrum by itself with some correlation factor Co(A;Ao). 



This can be deduced by considering the averaged form factor (4.9). The scaling function Cq should correspond to 
GUE statistics. Assuming that there are no cross-correlations, one obtains the result ( ^.5[ ). This result is consistent 
with GUE whenever G(fc) ~ 0. However, if G{k) — 1, then the correct GOE result is not obtained, implying that 
cross-correlations between POs of different winding-numbers are important. The identifications of the mechanism that 
leads to these cross-correlations is probably related to the geometry characteristics of diagrams which are considered 



25 



in the field-theoretical approach. For example, one may consider a PO that follows closely a loop-shaped orbit (say 
clockwise) and then a second loop-shaped orbit (say anti-clockwise) thus creating an 8 shape. However, due to the 
time reversal symmetry there may be a relate PO that follows the second loop clockwise rather than anti-clockwise. 
Thus, the two related POs will have different winding numbers, implying the conjectured cross-correlations. In this 
respect the GOE GUE crossover should be similar to symmetry breaking of spatial symmetry. However, the GOE 
— > GUE crossover cannot be described merely in terms of change of the scaling constant Ag Aq , since it involves 
also a change of the scaling function. Using different words, unlike spatial symmetries, the GOE — > GUE crossover is 
involved in changing universality class. 



B. Multi-cell billiards 



Complexity is associated with having non-trivial classification of POs. Disregarding symmetries, the simplest 
classification is induced by proximity. In the subsequent subsections we shall discuss the significance of spatial 
proximity for the analysis of multi-cell billiards. It should be emphasized that the main ideas can be generalized to 
cases where phase-space rather than spatial proximity is relevant. 

A multi-cell billiard is a structure that is composed of ballistic cavities (cells). See for example |2^,^. The cells 
are connected by small openings (holes). A one-cell billiard means a simple ballistic billiard. A two-cell billiard is 
constructed of two connected cavities. An iV-cell billiard with iV ^ oo is a prototype for a ID diffusive system. It is 
most convenient to assume that the last cell is connected to the first one, thus creating a closed chain. More generally 
one may discuss d-dimensional structures. Note that the dimensionality of the cells-structure may be lower from the 
actual dimensionality of space. For example, one may discuss a ID chain that is composed of 3D cells. For the time 
being we assume that there are no special spatial symmetries, neither reflection nor translational. 

For a simple ballistic billiard the only classical time scale is the ballistic one (average time between reflections). 
There is no distinct ergodic time scale. POs which are longer than a few reflections explore the whole volume. In case 
of a multi-cell billiard POs may be classifled by the volume which is explored. This classification implies that there is 
a distinct time scale terg for ergodicity. For POs whose period is t < tf.rg the explored volume is typically smaller than 
the total volume. For t^rg <C t all the POs are ergodic as for a simple ballistic billiard. The non-universal features 
that are associated with the above classification are manifest in the one-point statistical properties of the spectrum 
and we shall argue that they should also be manifest in its two points characteristics. 

First we discuss one-point properties of of the classical spectrum that corresponds to a multi-cell billiard. These 
are best described in terms of the so called classical probability to return Pci{t). By definition Pci{t)/n is the average 
probability density for finding the particle in the vicinity of its initial point, where Q is the physical volume. For 
simplicity we refer here to configuration space rather than to phase-space. Note that Pci{t) is related, essentially by a 
Laplace transform, to the Ruelle Zeta function. The trivial eigenvalue corresponds to the ergodic behavior Pci{t) — 1. 
For our multi-cell billiards there is a distinct time scale terg, such that for t < terg the classical probability to return 
is larger than 1. More specifically, from the definition it follows that Pci{t) « il/ild{t) where ^d{t) is the volume over 
which an initial distribution spreads after a time t. In particular Pci{t) = N initially {t ^ ^o): where N is the total 
number of cells, and to is the time it takes to escape into adjacent cells. The relation of the one-point properties 
of the classical spectrum to Pci{t) follow from the observation that the latter can be expressed as a sum over POs. 
Disregarding repetitions one obtains ||l^ 

KD{t)=tP,i{t) . (5.6) 



Note that the so-called Hannay and Ozorio de Almeida sum rule (2.3) is a special case that corresponds to the ergodic 
behavior. For t to this sum rule is approximately satisfied separately for each of the cells, and summation over the 
POs of all the cells indeed gives the correct result Pci(t) = N . 

We turn now to discuss two-point correlations of a complex system. The key observation is that POs may be 
classified by the explored volume. This claim is especially trivial iit -^Iq. Each PO explores one of the volumes that 
corresponds to one of the cells. The POs may be classified by the cell which is explored. There cannot be any cross 
correlations between POs that belong to different cells, since a deformation of one cell does not affect POs that belong 
to a different cell. The same argument can be extended to longer times {t < terg). For a given cluster of connected 
cells, corresponds a class of POs that explore all of these cells. Any other class of POs will have a statistically 
independent spectrum. It is sufficient that the corresponding clusters will differ by one cell in order to extend the 
above argument. Namely, deformation of this one cell will affect only one of the corresponding class-spectra, implying 
no cross-correlations. 

It follows from the above discussion that the form factor can be written as a sum of statistically independent 
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contributions, namely 



K{k,t) = J2 Kc.{k,t) = C^{k,t)K^{t) (5.7) 

a a 

Where Ka{t) is the (partial) diagonal sum over POs that belong to the a-class, Ka{k,t) is the contribution of these 
POs to the form factor, and by definition Ca{k,t) is the corresponding correlation factor. At this stage it is appro- 
priate to suggest a conjecture and a further working-hypothesis that constitute the basis for applying our formalism 
to the case of a multi-cell billiard. 



Ballistic-Like Correlations (BLC): 

The BLC conjecture - A class of POs that explore a volume Clg is char- 
acterized by non-trivial correlation scale which is A = Xo{t, Slg)- 

The BLC working hypothesis - A class of POs that explore ergodicly a 
given volume is characterized by non-trivial correlations with the same scaling 
function Co{k) as for a ballistic billiard. 



The BLC working hypothesis leads to the following explicit expression for the form factor 

K{k,t) = Co{kX{t,n„)) K„{t) (5.8) 

a 

Here fla is the volume which is explored by POs of the a-class. The BLC conjecture seems plausible due to the 
general way of thinking which has been introduced in section III. Accepting the claim that POs of a generic ballistic 
billiard are characterized by the universal correlation scale Xo{t,Q), it is most natural to claim that the total Q 
should be replaced by the explored O in the more general case of a complex system. On the other hand, the BLC 
working-hypothesis has no a-priori justification. Actually we would like to understand the circumstances and the way 
in which it fails. 

The strategy of the subsequent subsections is as follows. Subsections C and D demonstrate how the BLC working- 
hypothesis can be used in order to make a semiclassical computation of the spectral form factor beyond the diagonal 
approximation. This is done mainly for pedagogical reasons. The applications in subsections E and F serve as a real 
testing ground for our approach. We consider there a quasi- ID classically diffusive chain. In the limit N ^ oo the 
quantal form factor should correspond either to localized eigenstates or to band structures, depending on whether 
the chain is disordered or periodic. Our main purpose is to give an answer to the long-standing question, whether 
these results can be derived, at least in principle, from the SCTF. We demonstrate that the BLC conjecture by 
itself is valuable, since it gives the correct time scales. Furthermore, we try to make a quantitative theory, using the 
BLC working-hypothesis. The final result is not completely satisfactory. However, it demonstrates the feasibility of 
performing effectively such a computation in the future. 



C. One cell with a hole 



This subsection constitutes an introduction for later analysis of multi-cell billiards. The cell is a do-dimensional 

simple chaotic cavity (2 < do). Its volume is = £q". There is a small opening whose area is a^"^^. A generic 
orbit will escape eventually out of the cavity. Each chord of the orbit is roughly of length £o, and the probability 
to escape, once hitting the surface, is (i?o/ao)'^"^^- Therefore the average escape time is to = Lq — i}o/ao°~^. The 
particle is further characterized by its wavenumber k. Its velocity w = 1 is insignificant physically and it is used 
merely in order to translate units of length into units of time. Disregarding the small opening the Heisenberg time 
is t% = L^ PS Qok"^"'^- The dimensionless conductance go is defined in the usual way as the ratio of the Thouless 
energy 2'Kh/to to the level spacing 2TT?i/t^. Hence 

50 = 4Ao = {kaof'-' . (5.9) 

Thus go is simply related to the ratio of De-Broglie wavelength 2TT/k to the hole size oq. The De-Broglie wavelength 
is assumed to be much shorter than io as to allow semiclassical considerations. Actually, in order to have a valid 
semiclassical result, it should be smaller or at most equal to ao . 
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The classical probability to stay is found by solving the differential equation dp{t)/dt = —{l/to)p{t) yielding 
p{t) = exp{—t/to). The corresponding diagonal sum is 

ifi(i)=5]A,<5(t-L,) = i.p(t) . (5.10) 

j 

The above sum includes those POs that survive when the 'hole' is opened. 



D. Two cell billiard 



Here the classical probability to return is found by solving the coupled equations dpi{t) — — (l/io)(?';(^)~Pr(i)) Siiid 
dpr{t) — ''{l/to)(vr{t)~-pi{t)), where pi and Pr are the probabilities to find the particle in the left or in the right cells 
correspondingly ||50(| . Taking into account the normalization by the total volume, the classical probability to return 
equals 2pi{t) where the initial conditions are pi = 1 and pr = 0. Thus one obtains the result 

Pc/(t) = l + cxp(-2tAo) = l+p(t)' . (5.11) 

There are three classes of POs, those that explore the left cell (/), those that explore the right cell (r), and those 
that explore both cells (&). The partial diagonal sums should satisfy Ki{t) + Kr{t) + Kb{t) — tPci{t), and also 
Ki{t) = Kr{t) = tp{t). Hence we deduce that Kb{t) = t{l—p{t))^. The ergodic time for a two-cell system is simply 
terg — to- Non Universal features are expected whenever t is less or of the same order. 

At this stage we wish to test the BLC working hypothesis of subsection B. Using ( ^.8[ ) one obtains the following 
result for the form factor 

K{k,t) = 2xCo{kX{t,no))-tp{t) + Co{kX(t,2no)) ■t{l-p{t))'^ (5.12) 

We wish to discuss the validity of this result. 

Consider first the case of large conductivity 1 ^ .go- The non-trivial behavior is manifested in the classical regime 
where Co = 1 and K{k, t) — tPci{t). Disregarding the small time regime to the form factor looks like that of a ballistic 
billiard whose volume is 2Qq. However, one important feature is missed. The non-universal peak around to should 
be compensated later, else Rqmik,e —> 0) will be larger (compared with the universal result) implying that there is 
no complete repulsion for small energy differences. Thus we conclude that Co should be slightly different from the 
ballistic one. One may consider such an effect to be a consequence resurgence (see section HTF). 

We take now the other limit go -€1 1- Here the semiclassical formulation should not be trusted, since the De-Broglie 
wavelength is much larger than the hole size. Still, it is interesting to observe whether the qualitative behavior of 
the form factor makes sense. On the Heisenberg time scale the form factor looks like a union of two independent 
spectra that corresponds to the two cells. Note that such form factor implies that the repulsion on small energy 
scales is not complete. Due to the small hole that connects the two cells, there is a non trivial behavior of the form 
factor on the classical time scale to- This relatively large time scale corresponds to very small energy differences. 
Specifically there is a deep whose functional form is K{k, t) oc [1 — p{t){\—p{t))\. This deep implies a strong repulsion 
on small energy scale. Indeed, this is the expected behavior. However, quantitatively the result does not make sense. 
Firstly, the classical time scale to should be modified as in the case of tunneling. Secondly, the above formula predict 
over-repulsion (negative density for finding energy levels). 



E. Disordered Chain 



We consider a closed disordered chain that consist of N cells. The escap e tim e out of a cell will be denoted as 
before by to. The number of cells that are explored after time t is simply ^Jt/to- The ergodic time (known also as 
the Thouless time in the present context) is ierg = N'^to- We shall assume from now on that t <C t^rg, meaning that 
the finite size of the chain should not play any role in the dynamics. The classical probability to return is found by 
solving a diffusion equation. The detailed calculation is presented in appendix B. One obtains Pdit) = N/ y/2iT{t/to)- 
The POs may be classified by the volume which is explored. We denote by Kn{t) the partial diagonal sum that 
corresponds to POs that explore n specified cells. Looking on POs of length t, the probability for finding a PO that 
explores n cells should be proportional to Kn{t). It is obvious that this probability should depend on the scaled 
variable n/ \ftJTo. In appendix B we derive the scaling relation 




(5.13) 
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Above, /(r) is the scaled probability function. It should have the normalization 1. The function /(r) is determined 
analytically in appendix B, and it is plotted in Fig.^. However, the main features are quite obvious. It should be peaked 
at T = 1, with a Gaussian-like asymptotic profile f{x) ~ exp(— x^). For small x one should obtain f{x) ~ exp(— 1/a;^), 
which implies an exponential decay of K njt) as a function of t for small n. In particular for n=l it should be in 
accordance with the one-cell result (5.1C). 

At this stage we wish to test the BLC working hypothesis of subsection B. One obtains the following result for the 
form factor 



oo 

Kik,t) = NY,CoikXit,nnoj)-K4t) = nJ^^I 



to 



ngoto 



n=l 



27r 



t 



ngoto 



(5.14) 



Here Kq{t) stands for either the GUE or the GOE scaling function of RMT. Changing from summation to integration 
one obtain the scaling relation K{go,t) = NgotoK{t/t*) , where t* = gp^o and 



K{t) 



1 



dx f 



-K, 



1/2 



for 
for 
for 



r < 1, 
T < 1, 
1 < r. 



GUE 
GOE 



(5.15) 



This scaling function is plotted in Fig.^. The non-universal breaktime t* is neither related to the one-cell Heisenberg 
time t'jj = goto nor to the volume dependent Heisenberg time tn — Ngoto- Rather, it is determined by the disorder. 
The non-universal breaktime t* rather than the non-relevant universal time scale tn signify the crossover to a recurrent 
quasiperiodic behavior. Assuming that up to t* there is roughly a classical diffusive behavior, with diffusion coefficient 
V = il/to, one obtains the correct result ^ = goio for the so-called localization length. 




FIG. 6. Solid line - the scaled probability distribution of the explored volume for a disordered chain. Dashed line - the Gaussian 
distribution of the winding number for a folded periodic chain. 
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FIG. 7. The scaled form factor for a disordered infinite chain. Solid line - GUE result, Dashed line - GOE result. The thiner curves 
are obtained by employing the diagonal approximation, while the thicker curves are obtained by employing the BLC working hypothesis. 
The dotted line illustrates the correct asymptotic behavior. 

The present application demonstrates the success of the BLC conjecture, as well as the limitation of the BLC 
working hypothesis. The generalization of these ideas to localization in higher dimensions is discussed in [pT| . The 
BLC conjecture is sufficient in order to determine the non-universal breaktime t* and the corresponding localization 
length. Qualitative agreement with scaling theory of localization in one to three dimensions can be demonstrated. 
In particular, the present theory is capable of giving a qualitative semiclassical understanding for the existence of 
a critical regime in three dimensional systems. On the other hand, the quantitative analysis is unsatisfactory. The 
BLC working hypothesis overestimates the effect of (negative) correlations. We believe that the statistics is modified 
due to leaking of PCs, thus leaving holes in the classical spectrum. As a result there is some effective clustering that 
affects the functional form of the correlations. See also pi|. 



It is interesting to discuss the implications of translational invariance. Rather than having spectral statistics that 
correspond to localization, one should obtain spectral statistics that corresponds to a band structure. This problem 
has been discussed by |^,^. Here we wish to clarify the present limitations of the semiclassical approach, and to 
illuminate some confusion which may arise concerning the classification of POs. 

The one point properties of the classical spectrum are not changed when the disorder — > order route is followed. The 
relevant classification of POs is still by the explored volume. POs that explore the whole volume (n = N) constitute 
a negligible minority since we assume an effectively infinite chain for which t <C ierg- Still, in the case of a periodic 
chain one may further classify the POs in the class n = TV by their self-symmetry. Specifically, There are POs that 
are invariant to translations of = 1, 2... cells. It should be emphasize that as far as the unfolded chain is concerned 
these self-symmetric POs constitute a minority within a minority. On the other hand, once the classical spectrum of 
the folded chain is considered, the primitive versions of these self-symmetric POs become dominant in the classical 
spectrum. The class index in the latter case should be interpreted as the winding number. The distribution of the 
folded POs with respect to is a Gaussian p6| , p8| |. Let us denote by K^{t) the diagonal sum that corresponds to 
the class of POs whose winding number is v. Summing over u one should obtain tPci{t) where Pci(t) = 1 as for any 
simple ballistic billiard. Hence we deduce that 



where g{x) is a normalized Gaussian with unit dispersion (see Fig. ^). The distribution is peaked at I'—O, still the 
ly—Q POs become a minority with respect to 7^ if <o ^ ^■ 



F. Periodic chain 




(5.16) 
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From the above discussion it should be clear that the relevant classification of POs is totally different for the folded 
chain when compared to the unfolded chain. If one wishes to study the statistical properties of the band spectrum, 
one should consider an unfolded chain, and the corresponding classification is the same as for the disordered chain. 
The POs that determine the statistical properties of the spectrum for a given quasi-momentum, will not he involved in 
the semiclassical analysis of the band spectrum ! 

Thus, the only difference between disordered and periodic chain is associated with the two-point statistics of the 
classical spectrum. For a disordered chain the correlation scale is class-dependent, namely A — Xo{t,nno), while for a 
periodic chain A = Ao(t, f^o) for all the classes. Using the BLC approximation scheme, taking into account symmetry, 
one obtains the result 



K{go,t) = 




(5.17) 



The above analysis essentially coincides with that of [g8| for t < t^, where = goto is the one-cell Heisenberg time. 
Namely, for short times K{t) ^ \ft. Here we have tried to extend the semiclassical analysis beyond using the 
BLC working hypothesis. For t'jj < t we have obtained K{t) ^ 1/Vt, which does not agree with the correct ^ 1/t 
behavior. Thus, in order to have a proper semiclassical theory for band spectrum we should know the way in which 
the scaling function is modified. Again, we conclude that the BLC conjecture is effective in the determination of the 
various time scales involved, whereas the BLC working hypothesis is unsatisfactory. 



VI. NUMERICAL STUDY 



In this section we use our numerical database to study classical correlations in three different billiards: the 2D 
hyperbola billiard, the 2D Sinai billiard and the 3D Sinai billiard. One possible approach is to compute p{x) directly 
for various L-windows. A smoothing procedure is then essential in order to get rid of fluctuations. Our numerical 
strategy is different. We prefer to fix the L-window and to vary k. Our object of study is the correlation factor 
C{k,L) = Kci{k, L)/K]j{L). In order to avoid numerical artifacts we have not applied a smoothing procedure. 
Rather, the cumulative quantity is considered, namely 



^ik^ Lj) 1 2 



HKL)^ Tdfc'^^^ rdfc'%4^^ . (6.1) 

Jo 9eS Jo l^pll^q^pql 

Here J^p is over all the distinct lengths within a given L-window, and J^q over all the POs that correspond to 
a given length Lp. We have divided C{k,L) by g^Q which is the effective degeneracy within the specified window. 
This has been done in order to get rid of the trivial correlation factor which appears due to symmetries. This way 
the integrand always approaches 1 for k oo. Fixing the L-window and probing the spectrum within the specified 
window with larger and larger k is both mathematically appealing and numerically convenient. (It is less convenient 
to fix /c and to vary L). This is because our classical databases is limited to a narrow span of L values, while A: is a 
free parameter that can be set to any value. 



If the form factor complies with RMT predictions, then we expect (see subsection II F ) 



C(fc,L) « minfl,^) . J-nst/c-i , k<2n/XiL) 

' ' V Kd{L)J \ 1 , k> 27r/A(L) ^ ' 

and consequently 

const fc"* , k < 27r/A(L) 
const + k , k > 2tt/\{L) 

The above strictly holds for CUE, and holds approximately for GOE 
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A. The hyperbola billiard 



The desymmetrized hyperbola biUiard is a planar billiard that is defined by the curves a; = 0, y = 0, y = x, y = 1/x. 
The billiard was throughly studied by Sieber . The compliance of Kd with GOE for the desymmetrized hyperbola 
billiard was already demonstrated numerically in ||l9|. There a database of about 100000 POs with length L < 25 
(not including time reversal) has been used. Here we wish merely to demonstrate our numerical procedure using 
the same database |3^] and to illustrate the utility of I{k,L). The results are shown in figure ^, where we choose a 
few L intervals with 5L = \. The overall fc-dependence of the numerical data agrees quite well with the theoretical 
curves of GOE and is markedly different from the Poissonian (no correlations) case. This is a clear indication of the 
existence of classical PO-correlations that indeed conspire to give the correct GOE form factor. The agreement of 
numerical and theoretical data is also an indication to the correct L dependence, at least on average. The changes in 
L of the theoretical curves (GOE) between the lower and upper values of the L's considered are minute in comparison 
to the numerical dispersion, and therefore it is impossible to reach any definite conclusion concerning the detailed L 
dependence. 
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B. The 2D Sinai billiard 



Here we considered the 2D Sinai billiard with edge of length 2 and an inscribed circle of radius R ~ 0.5. The billiard 
was desymmetrized into quarter and POs were computed by Schanz and Smilansky [^8| . The database consisted of 
about 5000 POs in the interval L < 10 (not including time reversal and x y reflection symmetries). Our numerical 
results of /(fc, L) were computed in the interval 7 < L < 10. We found it very useful to separate the POs into classes 
according to the number n of their bounces with the circle. The actual justification for this separation into classes is 
a-posteriori, but following is also an a-priori motivation. Each bounce with the inscribed disc significantly increases 
(on average) the instability of the orbit. Hence, POs with different n's have in general very different amplitudes. 
Therefore, if we calculate C{k,L) with the amplitudes Aj taken into account, then the relatively few POs with 1-2 
bounces dominate and overwhelm the much smaller contributions of the other POs (which are much more numerous). 
On the other hand, if we ignore Aj^s (expect sign) and replace them by some average value, then orbits with small 
number of bounces will be effectively ignored. (This state of affairs is valid when L is relatively small and the n = 1 
class dominates. In general, the important factor is the multiplication of the density dn{L) of orbits with n bounces 
near length L with the typical amplitude A„{L).) In order to extract as much information as possible, it is imperative 
to separate the POs into classes according to n, such that all of the Aj in the same class are of similar magnitude. The 
numerical results are shown in figure |l^. Similarly to the case of the hyperbola billiard PO-correlations are clearly 
observed, and the general appearance is similar to the predictions of GOE that correspond to the length interval 
investigated. The numerical dispersion is, however, too large to allow a detailed test of the L-dependence. 




FIG. 10. The function I{k, L) for the 2D quarter Sinai billiard. POs were con- 
sidered in the length interval 7 < L < 10. 

The numerical results indeed indicate that the division into classes according to n is very useful as suggested above. 
It confirms the idea that there are non-trivial correlations between POs within the same class. This information is 
very valuable in the quest for the classical origin of PO-correlations. We shall exploit it in the case of the 3D Sinai 
billiard (see below). The existence of inter-class correlations as was shown above does not exclude the possibility of 
having intra-class correlations, but we did not investigate this issue. 

C. The 3D Sinai billiard 

In this subsection we investigate the classical PO-correlations of the 3D Sinai billiard. The 3D Sinai billiard is 
the natural extension of the familiar 2D Sinai billiard, and is defined as the free space between a cubic 3-torus and 
an inscribed sphere. It is presented in figure |l^. We consider a torus of edge 1 and sphere of radius R ~ 0.2. The 
billiard is desymmetrized into its fundamental domain (1/48 of the original billiard). The main reason for investigating 
this system is that we had an extensive database available that consisted of 7 • 10^ POs (not including time-reversal 
conjugates). This allows us to perform meaningful statistical checks. The 3D Sinai billiard is plagued with non-generic 
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neutral families of POs ("bouncing-balls") which were already mentioned in connection with the mixe d boundary 
conditions systems (subsection IV E ). As was discussed there, we should resort to the trace formula ( 4.33| ) that avoids 
the bouncing-balls, and takes into account only the generic, isolated and unstable POs. Consequently, the classical 
spectral density is defined with the coefficients Aj rather than with Aj . 




FIG. 11. The 3D Sinai billiard. Left: original billiard, middle: desymmetrized billiard, right: unfolded billiard. 



For a generic billiard system with time reversal symmetry one expects that the quantal form factor will satisfy 
Kqrn ~ Kqqe ~ 9eS^Dit) whcrc g^Q « 2. The similarity Kqm ~ Kqqe should hold globally while the similarity 
Kqm ~ g^gKi) holds for relatively short times. Also for the 3D Sinai billiard with mixed boundary conditions it 
has been observed |3|] that globally Kqm ~ Kcoe, provided that Kqm is normalized such that it approaches 1 for 
L oo. However, for short times Kqqe < Kqm < g^QKo- The deviation from GOE for short times is a non-generic 
feature of the system and will not concern us here. The deviation from the diagonal approximation will be discussed 
in the sequel in terms of classical correlations. 

As for the 2D case, we classify the POs according to the number of their bounces n from the inscribed sphere. The 
a-priori reason is that the amplitudes are very strongly correlated with n (much more than in the 2D case) and the 
division into the n-classes due to this hierarchy is plausible. In figure ^ we present the function /(fc, L) for n — 1, 2, 3 
and various lengths. The numerical results leave no doubt about the existence of classical PO-correlations. All of 
the data curves are strikingly different from the Poisson (lack of correlations) case. The division into classes turns 
out to be very useful in revealing very prominent correlations. Closer look at the data shows that the numerical 
correlations are generally weaker than GOE for low /c's (fc ^ 30) and are stronger otherwise. One might interpret 
this discrepancy as due to a non-generic effect (that prevails at least for low L's), but nevertheless we bear in mind 
that in the L-regime in which we worked the n = 1 term was dominant, and its correlations are not so different from 
the universal ones. Also, there is no clear and systematic L-dependence. The numerical data exhibit large dispersion 
that prevents us from concluding anything definite concerning the L-dependence. It is therefore useful to examine 
the correlations for a fixed L and various n's. This is done in figure where we fixed L = 5 and examined the 
correlations for n — This figure clearly indicates the systematic increase of correlations with increasing n. 
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FIG. 12. The integrated correlation factor 7(fe, V) for the 3D Sinai biUiard with R = 0.2 (desymmetrized) . POs are classified 
according to their number of bounces from the sphere. Top: n = 1, middle: n = 2, bottom: n = 3. The notation "L « 8" 
means a window in L near L = 8. On the left we present the data in linear scale, and on the right on log-log scale. The GOE 
predictions are truncated for k smaller than the smallest eigenvalue fci w 23.5. 
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To get a hint about the source of the PO-correlations, we have fixed L = 10 and n = 2 and varied the radius of the 
inscribed sphere. The resuhs are shown in figure We clearly observe the growth of the correlations as R decreases. 
This effect is qualitatively explained by the following argument. Given a PO of the full 3D Sinai billiard (on the 
torus), it can be described by listing the locations of the centers of the spheres from which it bounces. One should use 
for this purpose the unfolded representation of the billiard (see figure ^l]) . This list of "addresses" defines the code 
word W which is known to be unique [^l| , pT| . This uniqueness property does not hold for the desymmetrized Sinai 
billiard. Here, the code word W does not determine the PO uniquely. Periodicity is required in the desymmetrized 
domain, hence one has to determine additionally the symmetry element g according to which the PO closes on itself in 
the unfolded representation. This is illustrated in figure ^ For the 3D Sinai billiard, the symmetry group is Oh ]4^ ] 
which contains 48 elements. Therefore, there are up to 48 POs that visit the same spheres in the unfolded domain. 
We denote this set as ^{W). If the radius R is small enough, the lengths of all of the orbits in jiW) will be similar, 
and the orders of magnitude of their instability amplitudes are expected to be the same. The signs of the amplitudes 
within 7(W^) will differ, however, according to the parity of the symmetry element. Hence we expect a better and 
better cancelations within each 'y(W) as — > 0. This qualitatively explains the effect seen in figure 111. 
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FIG. 15. Two different POs witfi the same coding. The POs visit the same discs, but correspond to two different symmetry 
elements. The illustration is for the 2D Sinai billiard, but the same principle applies also to the 3D Sinai billiard. 



To test this argument we performed a numerical experiment, in which we gradually destroyed the possible correla- 
tions in a controlled way. The results are displayed in figure Two cases were considered: R — 0.2, L w 10, n = 2 
and R = 0.2, L « 5, n = 5. We first plotted the full correlations in which the coefficients AjBj are taken into account 



(see (4.33)). As in previous cases, strong classical PO-correlations are observed. We next disregarded the Bj factors, 
and re-ca culated lik, L) using only Aj^s. As clearly seen from the figure, the changes were minor. Next, we replaced 
the Aj^s with their signs ±1, disregarding their magnitude. Correlations were somewhat diminished, but one can 
say with full confidence that strong and significant correlations still prevail. This proves, that the coefficients are re- 
sponsible to part of the correlations, but more importantly, it proves that there are correlations between the actions, 
regardless of the coefficients. This gives a justification to our point of view taken in the theoretical part, in which 
we regarded the Aj^s as "slowly varying". The following step was to perform a controlled randomization of the ±1 
coefficients between ^{W) sets. We picked each set of orbits 7(VF) and either kept the original signs or reversed all 
of them. This means, that we kept the relative sign within each 7(11^). The definite result was that PO-correlations 
persisted! This proves, that the correlations between POs are largely due to correlations within the sets ^{W). To 
complete the argument, we randomized also the signs within each 7(VF), and indeed correlations were completely 
destroyed. The overall conclusion is that we were able to identify the grouping of orbits into sets with the same code 
word W but with different symmetry ^ as a prominent source of the classical PO-correlations in the 3D Sinai billiard. 
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The above investigation of PO-correlations was purely classical. It is therefore very interesting to enquire whether 
the quantal levels "fee l" the observed classical PO-correlations. To do this, we compared Kqm{k,L) (based on the 
left hand side of ( 4.33 )) to (7gg-ii'£)(L) (based on the right hand side of the same equation). The quantal from factor 
was computed from the exact quantal levels and velocities, and g^QK]j{L) was computed from the classical PO 
data (we did not assume any sum rule). The results are shown in figure |^ Note, that we integrated the resulting 
form factors in L to diminish fluctuations, and started the integration at L = 5 where the density of POs is already 
very high. Evidently, the classical diagonal approximation g^^Ku is consistently and significantly larger than Kqm 
in the {k, L) regime that we checked. Strong classical non-trivial correlations are therefore needed in order to "bend" 
g^QKjj down to the correct Kqm- Indeed, in the relevant {k,L) domain, the classical correlation factor C{k,L) for 
the dominant n — 1 contribution is w Q.Tg^^ (see figure ^2|). Quantitatively, if we multiply g^oKn by the numerical 
factor 0.7, then we get on the average compliance with the quantum results, as seen in figurejlT]. Therefore, we have 
demonstrated a case in which classical correlations are clearly manifest in the quantum mechanical spectral statistics. 
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FIG. 17. Comparison between the quantal and the classical form factors for 3D Sinai billiard with R = 0.2. 



Appendices 

A. EFFECT OF DEFORMATION ON PO-S STATISTICS 

Consider a billiard and a particular PO. Assume that the boundary of the billiard in the vicinity of one hitting 
point is deformed, so that the normal shift is b. Assuming no bifurcation, the length of the PO will become Lj 
Lj + 6j, where 5j = 2bcos{6j). The incidence angle 6j is measured with respect to the normal direction at the 
point where the PO hits the deformed surface. In order to prove this result let us define a coordinate x along the 
boundary (without loss of generality we consider a planar 2D system). The hitting points are (xi, X2, a;„). The PO 
constitutes an extremum of the length functional C{xi,X2, ■■■,Xn), thus SC/dxg — for the PO. For a deformation 
6C = {SC/5b)-b + '^^{S£/5xs)-{Sxs/Sb)-b = {5C/5b)-b. Hence, in order to derive the required result, it is legitimate 
to fix all the hitting points and to consider merely the length change of the chords due to the deformation. This final 
step is trivial. 

We now wish to consider the effect of general deformation on the two point statistics. Assuming that there are 
no bifurcations, the new lengths are Lj Lj + 5j. Treating the differences y = Sj—Si as realizations of inde- 
pendent random variables with probability distribution G{y), one obtains (Y^ A*AjS{x — (Lj—Li) — (Sj—Si))) = 
J dyG{y){'}2 ^i^jH^ ^ U ^ i^j^^i))) ■ The prime implies omission of the diagonal terms i = j. It is implicit that 
the sum is restricted to some L-window. From here one easily derive the result C{k) ^ [1 — G{k)-{1—C{k))], where 
G{k) is the FT ofG(y). 

The case of symmetry breaking due to deformation that preserves family-structure should be treated somewhat 
differently. Here it is convenient to write the change of length as Lj — > Lj + Sar, where a identifies families of 
the folded system and r distinguish between families of the unfolded system. Disregarding cross-correlations between 
families the double sum can be written as {J2a^ij^rs ^i^J^i^ ~ ~ {SarSas)))- The ij summation is 

over POs of the folded system that belong to the a family. This sum can be manipulated as before, note that it is 
essential to separate the r ^ s from the r = s terms, yielding J dy{J2ij J2rs^i^J^i^ ~ U ~ (Lj—Li)) • [SrsS{y) + 
{l—Srs)G{y)]) — J dy{J2ij -^i^jS{x — y — {Lj—Li)) ■ 9[5{v) + {9~^)G{y)]) ■ From here one easily derive the result 
C{k)^[l + {g-l)G{k)]-Co{k). 
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B. DIAGONAL SUMS FOR THE CHAIN PROBLEM 



We consider a chain that consists of N cells. The average escape time out of a cell is found as in section V-C, 
yielding 1/to ~ 2 x a^"^^ /flij. Here fig is the volume of a single cell, and Oq"^^ is the cross section of each of the two 
'holes' that connect the cell to the other cells. Let x be an integer coordinate that identifies the cells, and let p{x) be 
the probability for finding the particle in the cc-cell. This probability satisfies the equation 

^ ^ -^[2p(.)-p(.-l)-.(.+l)] « . (B.l) 

Above the difference equation is approximated by a diffusion equation. This approximation is valid on time scales 
such that p{x) is smeared over many cells. 

For an open chain of n cells {n < N) we solve this equation with the boundary conditions p{n/2) ~ p{—n/2) — 0, 
meaning that there is a zero probability to return once the particle escapes from the endpoints ('absorption'). We 
denote by p{x\x') the solution, where x' is the cell where the particle is initially prepared. In order to solve the 
diffusion equation one should define an appropriate set of orthonormal functions, which is in the present case: 

, _ \ ^/2cos{7rsx/n) for s = 1,3,5,... tu o\ 

' \ V2sin(7rs2;/n) for s = 2,4,6,... ^ ' ' 

Then we obtain the solution 

oo 

p{x\x') = ^e-(-'^'*)/(2„^*«V,(a;')^.(x) (B.3) 
The classical probability to return is simply 

f"/2 ^ / ^2,2 
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If the chain were closed then p{x\x') —tl/n for t oo, and consequently one would obtain the asymptotic ergodic 
result Pnit) 1. However, for an open chain the particle will escape eventually via the endpoints- holes, leading to the 
asymptotic result Pn{t) — > 0. Specifically, for very long times ri^to <C t one obtains Pn{t) ~ exp{—Tr'^t / (271^10)) . For 
intermediate times {t ^ n'^ta but still to <C i), one obtains Pn{t) ~ n/ ^2TT{t/to). The latter result is not sensitive to 
the endpoint conditions and is valid also for a closed chain. For short times {t of the order to) the above solution is valid 
qualitatively but not quantitatively. In particular, for n — 1 one should obtain the one-cell result pi{t) = exp(— t/tg). 

We turn now to the estimate of if„(t) which is the diagonal sum that corresponds to a class of POs that explore n 
specified cells. For a closed chain there are N classes that correspond to each n < N, and one class that correspond 
to n = iV. The sum Kjsfit) + J2n<N ^^ri{t) should be equal to tPci{t). For t < NHo we have the result Pci{t) = 
PN{t) — N / yj2TT{t/to), and furthermore the first term Kjsi{t) is negligible. For any finite t we may take the limit 
N —f 00 thus obtaining the sum-rule 

For an open chain that consists of n cells there are (n— n'+l) classes of POs that explore n' cells. Thus we obtain the 
sum rule 

n 

Y,{n-n'^l)Kn\t) = t-pn{t) . (B.6) 

n' = l 



The sum rule (pB.SD is actually a special case of (B.6), it is obtained by considering finite t, taking the limit n 00 
and neglecting {—n' + 1). Using the above sum rules we may obtain an explicit expressions for the functions Kn{t). 

At this stage it is useful to note the scaling properties of Kn{t). Looking on POs of length t, the probability for 
finding a PO that explores n cells should be proportional to Kn{t). It is obvious that this probability should depend 



on the scaled variable n/y't/to. Taking into account the sum- rule ( B.5 ), one obtains the scaling relation 



Knit) = ^] . (B.7) 



^2tt 
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Above, /(t) is the scaled probability function. It should have the normalization 1. Substitution of (Bj) into 



(B.6), and changing from summation over n to integration over x = nj ^/tjt^ one obtains {x—x')f{x')dx' 
^E^iexp(-^2g2/(2a;2)). Hence 



V2^^exp (-2^(^5)2 



s=l 



where, by Poisson summation 

F{x) = '^^pI^^^^")') ^ 2x^exp(-2xV) . (B.9) 



s— — oo 



The Function F{x) satisfies -F'(O) = ■\/7r/2 and F(oo) ~ and the corresponding derivatives are F'{fi) = —1 and 
F'(oo) — 0. It follows that Jp°° f{x)dx — 1 while xf{x)dx = \/Tr/2. The function /(x) is plotted in Fig.^. It is 
peaked around x ~ 1.2 and it vanishes for either a; — )■ or a; — > oo. 

ACKNOWLEDGMENTS 

We thank Michael Berry, Thomas Dittrich, Bruno Eckhardt, Shmuel Fishman, Martin Gutzwiller, Jon Keating, 
Holger Schanz and Martin Sieber for helpful discussions. The research reported here was supported in part by the 
Minerva Center for Nonlinear Physics of Complex systems, and by the Israeli Science Foundation. 



[1] 

[2] 

[3] 
[4] 

[5] 
[6] 
[7] 

[8] 
[9 

[lo; 
[11 

[12 
[13 
[14 
[15 
[16 

[17 
[18 

[19 



[20 

[21 
[22 
[23 



To obtain ( ^.(j ) one should use the SCTF with k replaced by k' , and multiply both sides by w{k'~k) where k is the center 
of the window. Then one should perform a Fourier transform in the variable k' to the variable L. 

The existence of a fimte bifurcation-free energy interval requires a further restriction over the periods Tj . This restriction 
can be realized by using an appropriate windowing procedure. 

H. P. Baltes and E. R. Hilf, Spectra of Finite Systems, Bibliographisches Institut, Mannheim (1996). 

J. P. Keating "The Riemann Zeta Function and Quantum Chaology", Proceedings of the International School of Physics 
"Enrico Fermi", Course CXIX, Ed. G. Casati, I. Guarneri and U. Smilansky (North Holland 1993). 
K.G.Andersson and R.B.Melrose, Inventiones math. 41, 197 (1997). 

M.C. Gutzwiller, Chaos in Classical and Quantum Mechanics (Springer- Verlag, New York, 1990). 
The notion of third entropy is explained in ^ p. 153 and p. 307. 
M. V. Berry and M. Tabor Proc. Roy. Soc. A349 101-123 (1976) 

Ya. G. Sinai in Ceometric Aspects of Functional Analysis Lecture Notes in Math. 1469, Springer, Berlin 41-59 (1990) 
P. M. Bleher ibid pp 60-89. and Duke Math. Jour. 74 45-93 (1994). 
p. Gaspard and D. Alonso, Phys. Rev. A 47, R3468 (1993). 

D. Alonso and P. Gaspard, CHAOS 3 (4), 601 (1993). 

J. H. Hannay and A. M. Ozorio de Almeida, J. Phys. A 17, 3429 (1984). 
T. Harayama and A. Shudo, J. Phys. A: Math. Gen. 25, 4595 (1992). 
M.A. Sepulveda, S. Tomsovic and E.J. Heller, Phys. Rev. Lett. 69, 402 (1992). 
H. Primack and U. Smilansky Phys. Rev. Lett. 76, 1615 (1996). 

R. Artuso, E. Aurell and P. Cvitanovic, Nonlinearity 3, 325 (1990), and references therein. 
M.V. Berry Proc. R. Soc. London A400, 229 (1985). 

E. Doron, U. Smilansky and T. Dittrich, Physica B179, 1 (1992). 

N. Argaman, U. Smilansky and Y. Imry Phys. Rev. B47, 4440 (1993). 

N. Argaman, F.M. Dittes, E. Doron, J. P. Keating, A.Y. Kitaev, M. Sieber and U. Smilansky, 
Phys. Rev. Lett. 71, 4326 (1993). 
S. Fishman and J.Keating, private communication. 



D. Cohen, "Periodic Orbits Breaktime and Localization", chao-dyn/9612014 (1996). 

D. Miller, "Toward Semiclassical Theory of Quantum Level Correlations of Generic Chaotic Systems" (1997). 
A.V. Andreev and B.L. Ahshuler, Phys. Rev. Lett. 75, 902 (1995). 
O. Agam, B.L. Altshuler and A.V. Andreev, Phys. Rev. Lett. 75, 4389 (1995). 
[24] B.A. Muzykantskii and D.E. Khmelnitskii, JEPT Lett. 62, 76 (1995); Pis'ma Zh. Eksp. Teor. Fiz. 62, 68 (1995). 



41 



E.B. Bogomolny and J. P. Keating, Phys. Rev. Lett. 77, 1472 (1996). 
R.Scharf and B.Sundaram, Phys. Rev. Lett. 76, 4907 (196). 
T. Dittrich, Phys. Rep. 271, 267 (1996). 

T. Dittrich, B. Mehlig, H. Schanz and U. Smilansky, Chaos Solitons & Fractals 8, 1205 (1997). 
B.D. Simons and B.L. Altshuler, Phys. Rev. Lett 70, 4063 (1993). 
N. Taniguchi and B.L. Altshuler, Phys. Rev. Lett 71, 4031 (1993). 
J.Keating, private communication. 
Z.Rudnik, private communication. 

M. Sieber, H. Primack, U. Smilansky, L Ussishkin and H. Schanz, J. Phys. A 28, 5041 (1995). 

M. Sieber, U. Smilansky, S. C. Creagh, and R. G. Littlejohn, J. Phys. A 26, 6217 (1993). 

H. Primack and U. Smilansky, Phys. Rev. Lett. 74, 4831 (1995). 

H. Schanz and U. Smilansky, Chaos Solitons and Fractals 5(7), 1289 (1995). 

M. Sieber, PhD thesis. University of Hamburg, DESY preprint 91-030 (1991). 

We thank Martin Sieber for periodic orbit data of the hyperbola billiard. 

We thank Holger Schanz for periodic orbits data of the 2D Sinai billiard. 

Harel Primack, PhD thesis. The Weiz mann Institute of S cience (1997, in preparation). 



P. Dahlqvist and R. Artuso, preprint, :hao-dyn/9601013 (1996). 
L. A. Bunimovich. Chaos 5, 349 (1995). 

M. Tinkham, Group Theory and Quantum Mechanics, McGraw-Hill Book Company (1964). 
P. Dahlqvist, J. Phys. A 28, 4733 (1995). 

T. Kottos and U. Smilansky "Quantum Chaos on Graphs" (preprint, 1997). 

E. Doron and U. Smilansky, Nonlinearity 5, 1055 (1992). 

F. B. Bogomolny, Nonlinearity 5, 805 (1992). 
M. Tabor, Physica D 6, 195 (1983). 

R. Bliimel and U. Smilansky, Phys. Rev. Lett. 64, 241 (1990). 

F. Haake, M. Kus, H.-J. Sommers, H. Schomerus, and K. Zyczkowski, J. Phys. A 29, 3641 (1996). 
U. Smilansky, S. Tomsovic and O. Bohigas, J. Phys. A 25, 3261 (1992). 
H. Schanz, PhD dissertation, Humboldt-Universitat zu Berlin (1996). 



42 



